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The work described in this final report is a continuation of the 
work started under NASA Grant 5288. The discrete approach to modeling 
vegetation originally developed under Grant 5288 has been extended and 
applied to a number of practical situations. 

Electromagnetic backscattering from a layer of vegetation over a 
flat lossy ground is studied. The vegetated region is modeled by dis- 
crete lossy dielectric scatterers, for which the dyadic scattering am- 
plitudes and orientation statistics are known. A method is developed 
to compute the backscattering coefficients from the vegetated layer. 

The technique is valid for scatterers having characteristic dimensions 
comparable to a wavelength. 

The problem is solved by finding the mean field in the vegetated 
region and then using it in conjunction with the distorted Born approx- 
imation (first order multiple scattering) to calculate the backscattering 
coefficients. The mean field due to a plane wave obliquely incident on 
the vegetation is obtained by finding an approximate solution to the 
Foldy-Twe.csky mean equation in the limit of small fractional volume. 

This approximate solution is obtained by employing a two variable ex- 
pansion procedure. An examination of the mean solution shows that from 
the viewpoint of the mean wave, the vegetated slab can be replaced by a 
deterministic anisotropic medium. The anisotropic medium becomes uni- 
axial when the scatterers are assumed to have orientation statistics 
that are independent of the azimuthal angle. 

The backscattering coefficients are then calculated by employing 
single scattering theory, in which the scatterers are assumed to be 

iv 


embedded in the equivalent anisotropic bulk medium. The procedure is 
valid when the albedo of individual scatterers is small, that is, when 
the scatterers are highly absorbing. Formulas for copolarized and cross 
polarized backs cattering coefficients are given. Numerical calculations 
for the backscattering coefficients as a function of incidence angle are 
presented. 

The material presented in the report has been published, in part, 
and has been presented at a number of technical meetings. Reference to 
this published material is contained at the end of the report. 


CHAPTER I 


INTRODUCTION 

The objective of this dissertation Is the Investigation of electro- 
magnetic backscatterlng from a layer of vegetation over a flat lossy homo- 
geneous ground. The vegetated region, or canopy, is modeled by discrete 
lossy dielectric scatterers that have prescribed orientation statistics. 

The motivation for this work has been the need to describe quantita- 
tively the effects of the natural environment on electromagnetic wave pro- 
pagation and scattering, and to relate radar return to the physical charac- 
teristics of the vegetation and the underlying ground. These characteristics 
can be used to determine the biomass and leaf area index of the vegetation 
and the moisture content of the ground. Information of this type is neces- 
sary as input data for crop yield models. 

Models have been developed to serve the above applications 
for vegetated terrain. These models have been constructed by replacing the 
vegetated region with a random medium whose statistical characteristics are 
related to the physical quantities of the medium. The random modeling 
techniques divide naturally into two types: continuous and discrete. In 

the continuous case, the random medium is modeled by assuming that its per- 
mittivity £('x) is a random process whose moments, such as the mean and cor- 
relation function, are known. In the discrete case, on the other hand, the 
medium is viewed as a collection of dielectric particles whose position and 
orientation statistics are given. In each, the medium statistics are used in 

conjunction with Maxwell's equations to calculate average quantities of 
physical Interest. 
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For the continuous case, exact equations for the Bean and correlation 
of the electric field can be obtained [Frisch* 1968]* These equations are 
known as the Dyson and Bethe-Salpeter equations respectively. Although almost 
Impossible to solve even under the most Ideal situations* they provide an 
exact formulation for the quantities .of Interest. Under appropriate physical 
circumstances* such' as media with small correlation lengths or small fluc- 
tuations In the permittivity* perturbation theory can be used to simplify 
both the Dyson and Bethe-Salpeter equations to a tractable form. In the 
case of the mean field* perturbation methods have been used by Keller [1962 
and 1964]* Tartarskll and Gertsenshteln [1963]* Keller and Karal [1966]* 
Rosenbaum [1971] and others. In active remote sensing applications these 
approximate mean field solutions have been used along with first order renor- 
malization or distorted Born approximation to obtain the backseat tering coef- 
ficients, [Rosenbaum and Bowles, 1974; Stogryn, 1974; Hevencr, 1976; Fung and 
Fung, 1977; Fung and Ulaby, 1978; Fung, 1979; and Zuniga* et al., 1979, and 
many others] . 

Another technique used to obtain the scattered field from a continuous 
random medium Is the radiative transport approach. Here the transport equa- 
tions are obtained in terms of the statistics of e(x) [Tsang and Kong, 1978]. 

In the case where the medium correlation length Is large compared to wavelength, 
de Wolf [1971], Ito and Adachl [1977] have developed multiple forward-single 

i 

backscatter techniques. 

In 1945 Foldy presented the first systematic probabilistic formula- 
tion of the multiple scattering of waves by collection of randomly distri- 
buted scatterers. In that paper he derived the mean field In a medium of 
scalar dipole scatterers. Lax [1951, 1952) generalized his treatment to 
resonant size scatterers and TWersky [1962, 1964, 1967, 1970* 1978] 


3 

employed the Foldy-Lax method to find the mean field and effective dielectric 
constant for arbitrarily shaped dielectric particles. Ishimaru [1978a] has 
obtained an approximate equation for the mean field in both the scalar and 
vector cases. Ishimaru [1978b] also has found solutions to the correlation 
equation by employing the diffusion approximation. 

The perturbation procedure applies when the fractional volume occupied 
by scatterers is small. The theory shows that the effective medium, as seen 
by the mean wave propagating thrcugh a collection of non-spherical particles, 
is anisotropic. Lang [1981], Lang and Sidhu [1983] used the Foldy-Lax method 
in conjunction with the distorted Born approximation to calculate the back- 
scattering coefficients from a slab of arbitrary shaped lossy dielectric 
scatterers. 

Another work, using discrete scatterers to model vegetation, has been 
done by Du and Peake [1969], and they employed single scattering (Born approx- 
imation) without introducing an equivalent medium. Thus, they did not take 
into account the decay of the incident wave in the vegetation. This limits 
their theory to a much lower frequency and thin layers of vegetation. They 
also did not take into account the underlying ground. 

Recently Tsang, et al, [1981], have used a vector radiative transport 
technique to analyze a slab of discrete scatterers. The results are similar 
to those Lang and Sidhu [1982] except that terms representing coherent wave 
effects are missing. The method is limited to low albedo particles where 
absorptive loss is the dominant mechanism. The two methods require different 
input quantities. For the continuous model the average permittivity of the 
medium is required in addition to the spatial correlation function of its 
fluctuations. In the discrete case, the scattering amplitude of the indi- 
vidual scatterers is required as well as the position and orientation statis- 
tics of the scatterers. The scattering amplitude can be obtained experimentally 
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or by electromagnetic modeling of the individual acatterer. The latter 
modeling procedure relatea the scattering amplitude to the physical di- 
mensions and to the dielectric properties of the acatterer. 

One of the advantages of the continuous modeling technique is that 
computed quantities of Interest such as the backacatterlng coefficient are 
obtained directly In terms of the average dielectric constant of the medium 
and the correlation function of its permittivity fluctuations. In the case 
of the discrete approach, the average permittivity, the correlation function 
of its fluctuations and backacatterlng coefficients are all obtained in terma 
of orientation averages over the scattering amplitudes of the individual 
scatterers. The scattering amplitude, in turn, is then related to the 
electrical and physical characteristics of the individual scatterer. Although, 
computationally, the discrete approach is more complex than the continuous 
method, the discrete approach has certain important advantages. One of these 
is that the average permittivity of the scattering medium is now a derived 
quantity rather than one which is empirically determined. Another dif- 
ference between the two approaches is that the continuous method does 
not permit cross polarized backscatter [see Tan, et al. 1980], to first 
order in the albedo, while the discrete theory predicts a first order 
contribution [Lang, 1981]. For these reasons, the discrete modeling 
technique may prove more successful in relating remote sensing signa- 

t 

* 

tures to actual physical characteristics of the medium. 

The methodology which has been employed in this work models vege- 
tation by dielectric discs (leaves) and uses discrete random media methods 
to calculate the scattering cross sections of interest. Lang [1981] assumed 
the canopy was thick enough so that ground reflection could be neglected. 



a * 
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The forest was thus modeled by a half space of dipole discs (400 MHz to 1 GHz 
regime). An equation for the mean field In the half space was derived and, 
from It, an equivalent dielectric constant for the leafy medium was obtained. 
Following this, the distorted Born approximation was used to calculate nhe 
backscatterlng cross sections of the leafy half space, Lang and Sldhu [1983] 
account for the effect of a flat ground by considering a layer of dipole discs 
over a lossy homogeneous half space. The distorted Born approximation was 
again used to find the backscatterlng coefficient of the layer of leaves. 

It was found that, the backscatterlng coefficient could be decomposed Into 
three terms: a contribution from direct backscatter; a return from a wave 

doubly reflected from the ground; and, finally, a direct-reflected component 
from a wave singly reflected from the ground. In this thesis, the method is 
extended to arbitrarily shaped scatterers having characteristic dimensions 
comparable to a wavelength (so-called resonance region) . The development 
is generalized by using a matrix formulation. 

This thesis has five chapters. Although the introduction to each 
chapter should provide the reader with an outline of its contents, 
we briefly summarize what will be done in each chapter. In Chapter II, 
Maxwell’s equations are recast in an operator form. The coherent field 
equation obtained by employing Foldy’s approximation will also be dis- 
cussed. The correlation of field will be found by employing the distorted 

i 

Born approximation. 

In Chapter III, the mean field and Green's function of a plane wave 
obliquely incident upon the vegetation is obtained by finding an approximate 
solution to the Foldy-Twersky mean equation in the limit of small fractional 
volume. This approximate solution is obtained by employing a two variable 
expansion procedure. An examination of the mean solution shows that the 
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vegetated slab can be replaced by a deterministic, anisotropic medium. The 
anisotropic medium becomes uniaxial when the acatterers are assumed to have 
orientation statistics that are Independent of the azimuthal angle. Also, 
the average complex dielectric constant of the medium Is determined from the 
average field in the medium. 

In Chapter IV, the model is extended to obtain the transverse spectral 
density of the field, and the resonant backscatterlng coefficients of the 
vegetated layer. Simple expressions for the copolarized and cross polarized 
backscatterlng coefficients are obtained In terms of the dyadic scattering 
amplitude of an individual scatterer. 

In Chapter V, a general discussion Is presented and numerical results 
obtained by modeling a forest canopy as a collection of lossy dielectric discs. 
The results are examined. 
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CHAPTER II 

GENERAL FORMULATION OF RESONANT BACKSCATTERING 

* 

- FROM VEGETATION 

In this chapter we present the general formulation for scattering 
from discrete random medium In the limit of small fractional volume. We 
adopt an approach which has been originally developed by Lang [1981]. 

First, we consider the problem of scattering of a time harmonic electro-* 
magnetic wave from N discrete identical lossy dielectric scatterers which 
have random position and orientation. Scatterers are considered to be in- 
dependent of one another, and as a result, neighboring particles are not 
necessarily aligned. Then we will consider single scatterers. In both cases 
we will obtain the necessary equations Jn terms of the transition operator. We 
will develop an approximate equation for the coherent field by employing the 
Foldy approximation [Foldy, 1945] . We will also obtain the macroscopic form 
of Maxwell's equations and the macroscopic permittivity operator which de- 
scribes the average behavior of the equivalent medium. 

Finally, the distorted Born approximation will be employed to obtain 
the correlation of the field in the equivalent anisotropic medium. , 

Problem Formulation 

We consider the problem of scattering of time harmonic electromagnetic 
waves from N discrete scatterers located in a volume V as is shown in fig. 1. 
The particles are identical and each has a volume V^, and a relative dielec- 
tric constant e f . We assumed that the relative dielectric constant of 
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background medium is e^(x). The only real restriction on the background 
medium la that It be constant inside V. 

The location of the 1th particle Is specified by the vector x^ f 
extending from an origin 0 to the center of that particle. The particle's 
center is located by the center of the smallest circumscribed sphere in 
which the particle can be placed. Although the particles are identical they 
have a rotation with respect to a fixed direction. The rotation for the 
ith particle is specified by (^■(G^,^) where and are polar and azimuth 
angles, respectively, with 0£© 1 £ff and 0<4» i <27r. 

The electric field £ and the magnetic field H obey 
Maxwell's equations 

-iu)e 0 e r (x)E - VxH - -J 


( 2 . 1 ) 


VxE_ - iu)jJ^H * -M 


where a time dependence e itllt has been assumed. In (2.1) J_ and M are, re- 
spectively, the vector electric current density and magnetic current density. 



Fig. 1. Distribution of particles within volume V 
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Nov let us assume that a particle located at the or9gln^^c^i^actlerlzed by 
the function U(x) , where 

1» * £ V 


U(x) p 

o. x<V p 


Using eq. (2.2), we express e (x) as 

JL 

N 

£_(x) " e a(*) + Mx) I 0(4 " *,* 04 ) 

1-1 


( 2 . 2 ) 


V x e R 


where 


and 


A(x) - e r (x) - e d (x) 


(2.3a) 


(2.3b) 


(2.3c) 


U(x,n) - U(R(fl).x) 

In eq. (2.3c) U(x,£) Is the function U(x) rotated by 13 , and R(fi) is 
a rotational dyadic. 

Equations (2.1) can be written In matrix form. 


(2.4) 


— m 

“iwe rt e (x)I -Vxl 

u r — — — 


E ' 


J(x)* 

Vxl 

- <m 

• 

H 


M (x) 

is the unit dyadic and 


* I» 

I*H - 

m 


Using (2.3a) in (2.4), we have for all x 



-iu)e_e . (x)I -Vxl ' 

0 d ~ - * 

N 

'-iaie 0 d(x)U(x-X 1 »i2 i )| 0 


V 



+ l 


0 



Vxl -1 u)Pq I 

1*1 

0 0 

. 


H 

• « 


J(x) 

M(x) 


(2.5) 


Equation (2.5) can be represented in the operator form [Felsen and 
Marcuvitz, 1973; Harrington, 1968] . 
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N 

(L- I v )¥(x) « g (2.6) 

1-1 x 


where L Is the operator descriptive of the field equations. The ¥ is a wave 
vector characterizing the field variables and the g is a wave vector de- 
scribing the excitation. 

Since L 1 b ai operator* it is necessary for uniqueness to state that 
¥ lies in a prescribed domain of the operator L, a remark that is equivalent 
to the statement of initial and boundary conditions on the elements of ¥. We 
list below the matrix form taken by the operator L and the wave vectors ¥, g 
and V for the electromagnetic field. 



-Vxl 

m 

Vxl 

- ±w ol 


(2.7) 


¥ 


E(x) 

H(x) 




1 w e n AlJ (x-X. » fl . ) I 
0 1 — i - 

0 


0 

0 

- 


( 2 . 8 ) 


(2.9) 


The matrix multiplication of equation (2.6) represents a dot product 
between the elements of the 2x2 matrix operators L, V and the wave vectdr ¥. 
The elements of L and are dyadlcs while the wave vector elements are 
vectors . 


At times it will be convenient to write 
¥(x) - ¥ Q (x) + V s (x) 

where 


(2.10a) 


sMSmWu . 



| 



GRKT'"' a « O 

OF PuOft QUALITY 





So<s>' 


\(X) 

m 


and f - 



V2>_ 


H/i> 


(2.10b) 


and VqQO is the solution to eq. (2.6) when no scatterers are present, i.e., 

Vt Q (x) “ g (2.11) 

* 

and ¥ Is the scattered field from the particles. 


Single Scattering Equation and Transition Operator 
Before proceeding with the N particle scattering problem, we will 
consider scattering from one particle located at the origin. Putting N«1 in 
eq. (2.6) with x^“0 and we obtain 

(L-VW - g (2.12a) 

where 



(2.12b) 


iwe 0 AU(x,0)I 0 


0 


(2.12c) 


The definition of L and g are given by (2.7) and (2.9) respectively. The 
eq. (2.12b) can be written as 


■ V'o + ^s 


(2.13a) 


where 



■So 



e 

— s 

*0 “ 


and 

ib m 
Y s 



^0 

* m 



h 
~ s 

• m 


(2.13b) 


and is the solution to eq. (2.6) when no scatterer is present, i.e.. 
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14» 0 - g (2.13c) 

We have used the $ notation instead of the Y for the field here to 
remind ua that there la only one acatterer present. 

If we use (2.13) in (2.12), we get 

- V\l> (2.14) 

* s 

From (2.14) we see that the term on the right, can be viewed as the source 
of the scattered field. We write 

g eq « V* (2.15) 

where g gq is an equivalent source term. Since we know that V«0 when it p V^, 
the sources g eq , exist only inside of the particle boundaries. 

It is more natural to think of the equivalent sources as being Induced 
by the incident field Because Maxwell's equations are linear, hence we 

an write 

8 eq * Tip 0 (2.16) 

The operator T in eq. (2.16) is known as the transition operator in the 
scattering literature [Lax, 1951]. Now using eq. (2.15) and eq. (2.16) in eq. 
(2.14) and multiplying through by l” 1 , we have 

ti> 8 - L”” 1 T \p Q (2.17) 

Thus the knowledge of T completely characterizes the scattering properties of 

1 

the particle and quantities of physical interest. 

The transition operator is a linear bounded, operator and, as a result, 
can be expressed in integral form 

g eq (x) - TiJ> 0 » j dx't(x,x'H 0 (x') (2.18) 


where the limits for the integral extend over all space and 
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t(x,x") 


hi -12 


-21 -22 


(2.19) 


As shown In Appendix A, for a dielectric scatterer (relative permeability - 1) 

a 

*21 " -12 " -22 " 0 


and 


( 2 . 20 ) 


in ' ^ 


The tte.aO In eq. (2.20) Is the same as the one Lang [1981] obtained 

-1 


when he solved the vector wave equation and 


iwu. 


is the. normalization constant 


associated with the source term. 

One can show that t Is 0 when x or x' are outside the particle 
[Frisch, 1968], i.e.. 


|(x»20 * 0 > x % V p or x' g V p (2.21) 

The property follows directly from the fact that the equivalent sources for 
the scattered field are located within the particle boundaries. 

The transition kernel t can be expressed in terms of its plane wave 
representation t, i.e.. 


t(x,x") 


— fdk dk't(k,k')e i( -‘- " K > 

(2ir V J 


or inverting eq. (2.22) 

t(k,lO 


(2tt) 3 


dx_ dx"t (x,x')e 


i (k * x - k' 




( 2 . 22 ) 


(2.23) 


In a similar manner, we can define the Fourier transform of t(x,x") 


t(k,k") 


(2tt ) ' 


dx dx't(x,x')e~ i( -*- * (2.24) 


We have used the notation that h is the Fourier transform of h. More 
specifically 
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» 

(k) “ dx h; 

i 


, , , -ik*x 
h(x)e 
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The dyadic scattering amplitude, f(£,i_), is defined in terms of the 


asymptotic expression for e g in the radiation zone. We have 


ik Q x 


i s <2L»i) * I (0,i) 


( 2 . 25 ) 


where k^ * w/y” Cq is the free space wave number, 

x 

0 is a unit vector in the x direction, 0 = - — r 

“ 1*1 

and i^ is a unit vector in the direction of the wave incident upon the 
scatterer. 

The relationship between f and t can be found for large |x| . (Appen- 


dix A) . The result is 


f(0,i) = 2ir (1-0 0) • |(k 0 0;k () i)*(I-ii) 


( 2 . 26 ) 


From this relationship, we see 


0-f = 0 


f*i = 0 


( 2 . 27 ) 


Thus f is a four component tensor — all combinations of two incident 
and two scattered polarizations. We also note that | does not completely 

determine t but rather only partially specifies it. 

Before concluding this section, transition operators for particles 
that are not located at the origin will be needed. As we know, the equiva- 
lent sources g^ for a particle located at can be related to the in- 
cident field. It follows that 


= T.<|- 0 = |t.(x,x0^ i) (xOdx' 


( 2 . 28 ) 


By shifting the sources and the incident field to the origin, ^ can 
be related to t. One finds that 
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t^x,*.') * tCx-X^x'-X^ (2.29) 

Note that throughout the discussion the dependence of t on rotations 
has been suppressed for convenience of notation. 

Approximate Mean Field Equation 

In this section, we will develop an approximate equation for the co- 
herent field by employing the Foldy approximation [Foldy, 1945]. The equation 
is in terms of the transition operator and the particle density of the 
equivalent medium. 

After the equation has been derived, it is pointed out that inside V, 
the mean field obeys Maxwell's equations with free space medium permeability 
and macroscopic permittivity that is inhomogeneous, anisotropic, and spatial- 
ly dispersive. 

The total field H 1 (x) can be thought of as a sum of the incident 
field fg plus a sum of the fields scattered from each particle, 

We have 

V = V n + l 4 ,U; (2.30) 

0 i=l S 

The total field incident on the ith particle is called th£ effective 

field and is denoted by Thus T^ ^ represents the equivalent sources 

generated by the incident field in the ith particle and thus the scattered 
field by the ith particle is: 

* L -1 T. (2.31) 

s I 

Using eq. (2.31) in eq. (2.30) we have 

N . ... 

Hx) = (x) + l L T. (2.32) 

o - i=1 i 

This is the equation that we wanted to obtain. Now, we average this equation 


and the result is 
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<Y> * y 


o + l L_1 <T i 

u i-1 1 


,Ci)> 


(2.33) 


To obtain an approximate equation for the mean, we follow [Foldy, 1945] and 
assume 

(2.34) 


^(i) ~ <y> 


This means that the random quantity f v ’ is to first order, equal to a de- 


( 1 ) 


terministic quantity. Using eq. (2.34) in eq. (2.33) and noting that 

<T ~ <x <4»>> * <T ><y> 

i i i 

We have the approximate equation for the mean fields 

N 


< 4 » 


—1 

> = ¥_ + 7 L <T.xy> 
0 i'l 


(2.35) 


(2.36) 


We need the dependence of average T^ upon and S2^. For this, we 
assume that the position vectors X. and rotation vectors are random 
variables that are specified by a 5N dimensional distribution function. 
Assuming the particles have identical distribution functions, we have: 


p x (x,«> = P x n (x,u) i = 1» 2, . . . ,N (2.37a) 

— i—i — — 

where w-(6,<J>). We assume that the particle’s location and rotation are 
independent, thus 


p x * p x ( 2 )p n ( “ ) 


with the usual property: 

P^(x)djx * 1 
V 


P Q Cu)du = 1 


The particle density is defined by 


(2.37b) 


(2.38) 


so that 



p(x)dx * N 
V 


(2.39) 


(2.40) 
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Now, we have 


<T ± > « <T(x l ,Sl 1 )> 


ds 


j j 
V 4tt 




where 


- J ds p x (s)T(s) 

V " “ 


T(s) 


j d “ P q(“) T (s,w) 

4tt 


(2.41) 


(2.42) 


In (2.42) the bar over T has been used to indicate only an average of angular 
variables. By substituting eq. (2.41) in eq. (2.36), noting that the scattered 
terms are identical, and by using eq. (2.39) we obtain 


<Y> = ¥ + 


ds p (s)L ^T(s)<¥> 


V 


Multiplying from the left by L and using eq. (2.11), we obtain 
$8<V> * g 


where 


S£ = L - 


ds_ P ( jO T ( s) 
V 


(2.43) 


(2.44) 


(2.45) 


This is the equation for the mean field which has been obtained essenti- 
ally by assuming that the effective field is approximately equal to the 
mean field (equation 2.34). Equation (2.34) is only valid when the frac- 
tional volume occupied by the particles is small compared to the total volume 
[Twersky, 1978]. We shall refer to a distribution of scatterers satisfying 
this condition as a sparse distribution. 

We now write the equation for the mean in a more explicit form. Using 
eqs. (2.28) and (2.29) in eqs. (2.44) and (2.45) we obtain the macroscopic form of 


Maxwell's equation: 
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L<H , (x)> - 


ds 


dX -P (8 ) t <X-S , X'-S ) < V/ (x - ) > ■ g 


where 


(2.46) 


t(x,x') * JdwP n (u)t(x,j{',jJi) (2.47) 

4ir 

Here the kernel t(x,x',w) is the same as eq. (2.29), however, we have explicitly 
shown its dependence on the angular coordinate w. 

By writing eq. (2.46) out in vector components, we explicitly obtain 
Maxwell's equations. 

The macroscopic form of Faraday's law: 


Vx < E^(x) > 35 1 wPq<H(x)> - <M(x)> (2.48) 

and the macroscopic form of Ampere's law with displacement current: 

Vx<H(x)> * “iu)<D (x) > ■+ <J_ (x) > , <D> = E-E <E> (2,49) 

U"“tJ 

where e e is a macroscopic permittivity operator which describes the average 
behavior of the medium. We have 


e (x) 

asp ’ — ' 




If J 

K o V 


dx'p (£)| (x-^x'-s) * 


(2.50) 


This expression simplifies to Ie^x) when x y. V. To see this we note that 
when x £ V, we have x~s. f. V , since s_ eV except for a small region near the 
boundary. Now using eq. (2.21) we have t*0. When xeV, eq. (2.44) does not 
simplify in general. It describes anisotropic, inhomogeneous, spatially 
dispersive medium. 

Lang [1981] showed that eq. (2.50) reduces to some more familiar ex- 
pressions in certain special situations. 
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The C o rrelation Field Equation 

In this section we will calculate the correlation of the field. The 

distorted Born approximation will be employed. It is a single scattering 

approximation where the scatterers are assumed to be embedded In the equiva- 

* 

lent anisotropic bulk medium. The procedure is valid when the albedo of a 
single particle is small. The latter condition implies that the energy 
absorbed by a particle must be much larger than the energy scattered by it. 

We start by considering a volume V of equivalent medium which was 
mentioned before. There are N particles embedded in V as shown in fig. 1. 

The scattered field due to the ith particle can be calculated by modifying eq. 
(2.31). We assume that the incident field on the ith particle is the mean 
field <4*> and that the operator L is replaced by the equivalent medium 
operator se as given in eq. (2.45). 

We have 


Y 

s 


N 

l 

i-1 


v (i) 

se 


where 


(2.51) 


= ^“ 1 T,<4'> 

se i 

Before proceeding, we point out that our main interest in finding the cor- 
relation of the field is to use it to calculate backscattering cross 
sections. We define the correlation of the field fluctuations as 

Y- - Y - <Y > <4* > = 0 (2.52) 

f s s f 

Now computing the correlation of the fluctuation component of the scattered 
field, we obtain 


<4'^ (x)Tj(x)> = < 4' g (x)4'g (x)> - <4' s (x)><Y^(x)> 


(2.53) 


where 4 1 denotes the conjugate transpose of 4'. 
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4 » 

In eq. (2.53) i|> Is conformable to Y for matrix multiplication. The 

4 * 

elements of ¥ and ^ are vectors. Therefore, eq. (2.53) represents a 2x2 
matrix with elements of dyadics. 

+ 

Putting eq. (2.51) in eq. (2.53) and noting that a portion of <! f > 

o 8 

cancels the term <¥ ><V >. We also assume that N»l. Then, we v'btain 

s s 


<Y f (x)¥^(x)> - Jdwp n (w) < 'i , f (x)f|(x)> w 

47 ! “■ 


(2.54) 


where 

<V f (x)’fj(x)> w * | dnp (£)^ s (x.jOYg (x,sj (2.55) 

“ V e e 

with 

^ (x.s) ‘ 1 T(s}*¥> (2.56) 

e 

Here we have separated the average into rotation and coordinate space 
averages, and thus we introduce the conditional expectation, » wit h 

respect to given ok 

The field S' defined in eq. (2.56) represents the field scattered 
s 

e 

by a single particle located in the equivalent medium at s _. V is the volume 

of the medium. We see that equation (2.55) is just an incoherent addition of 

the single scattering contributions from each particle. This is a result of 

the assumed independence of particle statistics. 

The correlation of the fluctuation field ^(x) is calculated instead 

of the correlation of S' (x) , so that coherent effects will be eliminated from 

s 

the backscatter coefficient expressions. 

The backscattering coefficients are then directly related to the trans 
verse Fourier transform of equation (2.55) with respect to x and x^. evaluated 
at the upper interface (z=0). This work will be done in chapter IV. 
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To write equation (2.56) more explicitly, we introduce the Green v s 
function 0 (x^* jc") for the operator It satisfies 

J^G(x,x') * 1 6 ( 2L“"iL”* ) (2.57) 

where 2 is given by eq. (2.45). The I is a square diagonal matrix of the 
unit dyadics, 

. r 

I 0 


X * 


(2.58) 


and 


G(x,x') 


OilC*** - *) £ 10 (x,:>O 


12 


£ 21 ( x »>0 (*»*'’) 
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(2.59) 


where G.^, -i'" £ 2 1 an< * =22 are Green’s functions. [See Felsen and 

Marcuvitz, 1973]. Now equation (2.56) becomes 


¥ fl (x.s) ■ 
e 


d^G(x,5c") 


dx^t(x‘'-s,x"-s)<Y (x")> (2.60) 


Equation (2.60) simplifies in the low frequency or Rayleigh 
limit. 


Before proceeding to a specific application, we would like to point 
out the relationship between the distorted Born approximation as presented 
here and the Twersky equation for the correlation. Basically, the vector 
analog of the scalar Twersky equation [Twersky, 1964] can be found by applying 
the method of smoothing when the fractional volume is small. If one then 
solves this equation under the assumption of small albedo, the result should 
be the same as the distorted Born method we have employed. 


CHAPTER III 


MEAN FIELD AND GREEN’S FUNCTION IN THREE MEDIA 

*» 

The mean field and Green's function in three layered media are develop- 
ed in this chapter. We begin first by obtaining the transverse wave vector, 
then we obtain a fluctuation equation. Using the two-variable expansion procedure, 
the equation will be solved approximately. The two-variable expansion procedure 
is a perturbation technique for approximating solutions to differential equa- 
tions, and it is the principal mathematical technique used in this study. 

Some particular applications of this method has been made by Tsang and Kong 
[1976], Tan and Fung [1979], Tsang and Kong [1979], and Zuniga and Kong [1981], 

We conclude the mean field study by obtaining the mean field of three 
media in a small volume limit (pV^cd) . The small volume limit or sparse 
distribution is natural for vegetation where leaves and stems make up a re- 
latively small percentage of the volume. This in turn leads to the concept 
that the mean field acts as though it propagates in an equivalent medium 
which is anisotropic. Because of this condition, we expect interface effects 
to be relatively unimportant. The most important effect is the decay of the 
mean field as it progresses through and reflects from the ground. 

Finally, we will compute the transverse dyadic Green's function in 
three layered media but will only give explicit expressions for it in the 
slab region. 

The Evaluation of the Eigenvectors 
A. Transverse Fourier Transform 

To illustrate the application of the method developed in the previous 
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sections, we will calculate the mean field in a three layered medium. As 
shown in fig. 2, this consists of free space at the top, the vegetation layer 
and the ground. A plane wave is assumed to be incident at an angle 6^ and 
scattered at an angle 8 g . The unit polarization vectors for incident and 
scattered wave are shown in fig. 2. The interface between the ground and the 
vegetation is taken to be smooth. 



Fig. 2. Geometrical configuration of the problem 


Medium G is free space having permeability and permittivity e^. 
The scattering particles are all assumed to be identical and each has a 
volume Vp, relative dielectric constant e r , and free space permeability p^ in 
Medium 1. Medium 2 is the ground underlying the vegetation with relative 
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dielectric constant e and permeability p . 

S ^ 

The equations satisfied by the mean fields in the three media, which 
are mentioned abo^e, have been derived in chapter II. There it was shown 
that the average electric field, <jS(x)>, and the average magnetic field, 
<_H(x)>, obey Maxwell’s equations in an equivalent medium having relative 
permeability 1 and -relative permittivity e (x). The relative permittivity 
e (x) is dyadic since the scatterers are arbitrarily shaped. 

In this section we consider first the abstract formulation of the 
guided-mode wavevector representation of Maxwell’s equations. We start by 
writing the macroscopic form oi Maxwell's equations which were obtained in 
chapter II. Assuming g=0, we have 


L<y(x)> - 


ds 


dx'p (s) t (x-s , x^-s.) <¥ (xO > 


(3.1) 


where 



H* 

" iwG O e d(~)i 

-Vxl" 

ac 

| 

L = 



with e (x) = e,(z) * < 


Vxl 

1 

Mil 

O 

3- 

3 

*H 

1 

( 


1 , z>-d 


e , z<-d 
g 


and 


<V(x)> - 


<JE (x) > 
<H(x)> 


t = - 1 


iuip 


0 


(3.2) 


(3.3) 


(3.4) 


0 0 

The left side of eq. (3.1) becomes zero when x/V . To see this we note that 
when x^V, we have x-s>V p since seV. Now using eq. (2.21) , we have t=0. Thus the 
integral over V in eq, (3.1), can be replaced by infinite limit if x is 
restricted to lie in the slab region. 


L<4'(x)> 


| | ds| dx' 


p (s) c (2L~ s.* ”*s_) <Y(x')> > rectj(z) 


npr'*^2f> 3 n rr> 

*** ' '■ ■; r 
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(3.5) 


where 


0 , z>0 


(free space) 


rect^(z) - J 1 , -d<z<0 (slab) 


0 , z<-d 


(ground) 


(3.6) 


The general field in the vegetation can be described either in terms of a first- 
order system of field equations, or it can be reduced to a higher-order equation in 
terns of a particular field component. The reduced field formulation (eq. (3.5)) 
frequently leads to analytical complexities in identification of energy express- 
ions, reciprocity properties, eigen modes, etc.. The first-order formulation 
avoids many of these difficulties. In this work the first-order formulation will 
be used. Also this formulation will be more suitable to apply to the two-variable 
expansion procedure. 

In order to simplify eq. (3.5), we should like to take the transverse 
Fourier transform of it. For this purpose, we assume p has no transverse 
variation, i.e., 


P (s ) - p(s) 


(3.7) 


where 


s * s + sz l 


Also we define 


X(x t ,x^,z-s,z'-s) « d^tO^-s^x^-s^z-s.z'-s) 


(3.8a) 


and let 


-t " — t ~ -t 


ds_ = ds. 
— t — t 


Using eq. (3.8b) in (3.8a) we have 


(3.8b) 


X (2^ »2Lt ;z ~ s » z ds^t (-s t ,xj>x t -s t ,z-s ,z '-s) 




(3.8c) 
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As seen in eq. (3.8c) we suppress the dependence by using eq. (3.8a). 
Substituting (3.8c), and (3.7) into (3.5) we have 


L<4'(x)> 


in 


dj*"p (s ) X (x^-xT ; z-s , z "-s ) <4 J (it ) > 


rect ( j(z) 


(3.9) 


In eq. (3.9) we have written the difference in transverse coordinates because 
the particle density p does not depend on £ t> Because of the exponential 
dependence of the incident wave on x and the invariance of the mean equation 
in the transverse direction, we assume 

ik f * x*. 

<1'(x)> = <'i'(Ji c .z)> * <4'(k t ,z)>e (3.10) 

where 

k t = k^yP , k t = kgSin0 (3.11) 


The system of ordinary differential equations is obtained by substituting eq. 
(3.10) into eq. (3.9). This is equivalent to taking the transverse Fourier 
transform of eq. (3.9). We find: 

,z ')>i rect d (z) (3.12) 


L<4'(k t ,z)> 


I 


ds |dz'p (s ) x (kj. , z-s , z ■'-s ) <¥ (k^ 


where 


<Y (k. ,z)> * 


cE(k t ,z)> 


<H(k t ,z)> 

The transverse Fourier transform of the operator L is 


L = -iK + r-f- 

dz 


with the component operators K and T defined by 



We 0 e d (z) * 

txl’ 



0 

«• 

0 T 
~z_ xl 

K * 

* 


and 

r = 




-k^xl 

“V. 



z°xl 

0 

m 


(3.13) 


(3.14) 


(3.15) 


From eqs. (2.20), (3.8a) and (3.11) we see that 


f (k t ,z-s,z'-s) r f(z-s,z'-s) 


TuCz'S^^-s) 0 


(3.16) 


B. Eigen-Functions for the Particle Free Medium 

The complete solution to eq. (3.12) Is obtained by a perturbation tech- 
nique under the assumption that the fractional volume is small. As a first 
step in implementing this procedure the homogeneous solutions to L<Y>*0 are 
found. The homogeneous solutions are of the form: 


<Y> = Y e 

a 


Ik z 
a 


(3.17) 


where Y and will depend upon whether z>-d or z<-d. Substituting eq, 
(3.17) into the homogeneous equation L<Y>*0 yields 


LY - (K-k T)Y - 0 

a a a 

or 


KY » k TY 
a a a 


(3.18) 


where K and F are defined by eq. (3.15). It is seen that the k and Y are 

U U 

eigen values and eigen vectors respectively of the operator K. 

To establish the orthogonality of the eigen vectors, we introduce the 

product 


(Y ,Y a ) - Y Y * E “E_ + H *H_ 
' a* 8 a p -a —8 —a —8 


(3.19) 


T + 

where Y is the transpose of Y. The adjoint operator L is defined by 


< L VV * <V L V 


For the product defined in eq. (3.19), we conclude 


<L V y g ) - (L v s . »+ L 1 L + » 8 
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Nov the adjoint eigen function problem for eq. (3.18) can be defined as 
4* + *4" 

k t a - K a r * a 

+ T + T 

Since K * K = K and T ■ T * T the adjoint problem is the same as the original 

+ + 

eigen value problem; thus ^ a m ^ a and K a m < a ' As a result we say the eigen value 
problem is self adjoint. 

We can now show that two eigen vectors having different eigen functions 

are orthogonal. Consider eigen vectors V a and and their corresponding eigen 

values A^ and A 0 . We have 
a p 

< K W - (K W ‘ 0 

«vv - <w - 0 

<WV - <WV * 0 
V r W - WV - 0 


V <r VV 


If A then 
a f 0 


< r W a 0 


(3.20) 


The eigenvectors^ and are orthogonal to each other as defined by eq. (3.20). 

Actual computation of the eigen vectors and values proceeds in a direct 
manner from eq. (3.18). The calculation follows the method of Felsen and 
Marcuvitz, 1973 and thus will not be repeated here. The results of the calcula- 
tion shows that there are four eigen vectors and four eigen values. The eigen 
values are given by 

k + (z) ■ 


±<(z) 


(3.21a) 
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where 


k(se) ■ 


K 

K 


z>-d 

z<-d 


with 


k * k^cosS , 


A. 2 - k 2 


k = k n /r 

g 0 g 


(3.21b) 


(3.22) 


g g t 

and k fc is given in ‘eq, (3.11) . We note that there are only two distinct eigen 
values; each has multiplicity two. Thus two eigen vectors correspond to each 
distinct eigen value. The eigen vectors are: 


'l' +1 (z) IB 


r4<‘> l 

L-n(z)h° J 


^±2 ( z) 


-h° i 
. n(z)v + (z)J 


(3.23) 


where the unit polarization vectors are given by 
. 0 _ 


0 , V 

v, (z) 


, u u 

V* 


(3.24) 

0 



r± 

2>-d 


0 


(3.25a) 

-g± 

z<-d 



with 


v ° ± - (+kx° + k^ 0 ) /k Q 


and 


0 


n(z) 


(+K + k.z^)/k 

g“ tr- g 


I Vd (z) 

V ^0 


(3.25b) 

(3.25c) 

(3.26) 


An examination of the eigen vectors shows that they obey the following 
orthogonality relationship 

(VV - 2N a 6 ug » a ’ e £ {±1 ’ ±2} (3.27a) 

where <5 0 is the Kronecker delta function and N is a normalization coefficient 
a3 a 


given by 
N 


± i ' iV 2) 


R.U) - 

1 Kq 


1 - 1,2 


(3.27b) 
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This result embodies the orthogonality relationship given in eq. (3.20), 
plus the fact the eigen functions corresponding to the same eigen value were 
choosen to be orthogonal. 

The eqs. (3.23) represent the eigen vectors for z>-d and fqr z<~d. For 
further use we write these eigen vectors explicitly in each region. We have 



for z>-d 


for z<-d 


(3.28a) 


(3.28b) 


with n = n /T . 

8 8 

In obtaining eq. (3.28b), we have used v^=v^ . Due to radiation con- 

“8 “ 8 - 

dition, v^ will not come into our development . 


The Solution of the Mean Wave Equation 
A. Reduction to Slowly Varying Coefficients: 

We start with substituting eqs. (3.23) in eq. (3.17), then we obtain 

<nz)> = (V_ ± a' + V_ 2 a")e~ iKZ + (H^b" + V +2 b")e' Uz (3.29) 

where the form of the eigen values and vectors will depend on whether z<-d 
or z>~d. 


The problem is now reformulated in terms of slowly varying parameters 
[Kohler and Papanicolaou ,1973], These parameters are defined by the trans- 
formation 


A(z) = a'(z)v^ (z) + a" (z)b? 
B(z) = b"(z)v° + (z) + b~(z)h° 


(3.30) 


Here A(z) and B(z) are the slowly varying wave amplitudes associated incident 
and scattered wave, respectively. 
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y (z) * (z) + f_ 2 ^° 

(3.31) 

y + (z) « » +1 v J (z) 4- v +2 h° 

where and V +2 are defined by eq. (3,23). The y’s are 2x1 matricies 
whose elements are vectors, therefore y” and y + are 2x1 matricies where 
each element is a dyadic. 


Using eq. (3.30) and (3.31) in (3.29) we obtain 

<V(z)> - y~Ae~^ K£ + y + Be lKZ 

Defining the vector $(z) by 

A(z) 


$(z) 


B(z) 


and the 2x2 matrix y(z) of dyadics by 

y(z) * [y (z)e 1KZ ,y + (z)e iKZ ] 

We cast eq. (3.32) in matrix form as follows: 

<i(z)> = y(z)<}>(z) 

where we denote the matrix y(z) in the two regions by 

y , z>-d 


y (z) = 


z<-d 


(3.32) 


(3.33) 


(3.34) 


(3.35) 


The equation (3.35) represents the transformation from field quantities 
to slowly varying amplitudes which we wanted to obtain. 

Using eq. (3.23) in eq. (3.31) we obtain 
l"(z) 


y (z) 


n (z)i°_xi (z) 


y + (z) = 


X + (Z) 


0 -f 

n(z)i ,xl (z) 


(3.36) 


where ,i + (z) are unit vectors in the direction of the up and down going waves. 


They are defined by 


£(z) - 


CiSt ± K £ )/ k 0 


(kt ± )/k g 


cassm. ; rs 
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(3.37a) 


I”(z) * Ih^h^ + v^(z)v^(z) 


(3.37b) 


Here I“(z) have the property l'"(z)'A(z) ■ A(z) and I (z)*B(z) ■ B(z). 

The I~(z) and I + (z) are unit dyadic tensors associated with up and 
down going waves respectively. 


Using eq. (3.36) in eq. (3.34) we obtain the y matrix as 


y(z) = 


l (z)e 


n(z)((Ar(z))e _iKZ 


I + (z)e lKZ 


n(z) ((i^xI + (z))e iKZ 


(3.38) 


One can define dyadic inner product similar to eq. (3.19) [Friedman, 1962], as 


(y + »ry + ) 


+ T + 
(y ) r(y-) 


(3.39) 


Substituting eqs. (3.36) and (3.15) into eq. (3.19) and proceeding with 


the matrix multiplication, we obtain the results 


+ T + + 

(y“) r(y-) = + 2N_I’ T ’ 


(3.40) 


+ T , +, 


with (y'~) T(y ) = 0. The normalization coefficient N’ T is given by eq. 


(3.27b). 


Now we define the inverse matrix, y , as: 


-1 , V r -1 - , v iKZ ,1 "K \ -iKZ.T 

7 c*> ■ y (z * e * + 2 n~ y (z)e i 


(3.41) 
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-1, . 1 
y (») * 25 


-l’"(z)e 


T| _+. . ~i<z 

I I (z)e 


-n(z)(i°xr(2))e iKZ 


/ >„0 _+/ s V “iKZ 
n(z) (i|.xl (z))e 


The inverse matrix has the property that 


y -1 ry - i. 


where 


I 0 

t= 


♦ + 

0 I 


Here I. is a identity matrix in the a space where I <{>“ 

Next we substitute eq. (3.35) into the transverse mean wave eq. 


(3.12) . We obtain 


Ly(z)<}>(z) 


^ jdsjdz" 


p(s)x(z-s,z"-s)y(z 


-)4>(o) 


where 


Ly(z)nz) - [(-iK+r~-)y(z)]^(2)-(-iKyfr^y)*+ry^- (3.4< 

Using eq. (3.34) we obtain 

(-iK+r~)y * (-iK+r~) (y“e~ iKZ , y + e +iKZ ) 

= [(-iK-iKr)y‘e" iKZ , (-iK+ir K )y + e +iKZ ] (3.4; 

Using eq. (3.31) in eq. (3.47) and proceeding the matrix multiplication, we 
obtain 


(-iK+r-j^-)y(z) = 0 


Hence eq. (3.45) becomes 


dd> (z> 


= ^ JdsJdz'p(s)x(z-s,z'-s)y(z')4>(z') j 


rect ( j(z) 


(3.49) 


, onmmi pags is 
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Multiplying both aides of eq. (3.49) by y *(z) and using (3.43) we 
find the fluctuation equation to be 

~tfr(z) -^jds|d*^p(s)y“ 1 (z)xy(*0*(*^)^ r«ct d (r) (3.50) 


with boundary conditions: 1) tangential components continuous, l.e., ly(z)$(z) Is 
continuous at z ■ 0* and z">-d, and 11) radiation condition at£«. Ve also have 


y“ 1 (*)fy(z") p^(*,*A«) - 


with 


-♦ -♦ 1 

v<+»-) f<+»+) 


(3.51) 


y(o,B)y + .\ B ( 5 ) t®.v (it ,t^ “iK((U“6z ) 

x u»* »*>■ ^ i -u^ 0 * * z *' i * » 

o,6 e {+»-) 


B. Application of the Two-variable Procedure: 

The eq. (3*50) can be solved along with the appropriate boundary con- 
ditions by using exact or approximate methods. The two-variable perturbation 
procedure will allow us to determine the approximate solution to eq. (3.50) 
directly without having to find the more complicated exact solution first. 

The main advantage of the two-variable procedure is the simplicity of 
the formalism. The higher approximations are more easily calculated. A 
disadvantage of the method, however, is that the proper choice of fast and 

slow variables are not always obvious [Cole, 1968]. 

We assume a slow variable z"6z, and then we assume p(z)»5?f(z). Thus 
pTz) Is a slowly varying function. Next we consider $(z) to be both a 
function of z and zT. We have 

00 

$(z) 2 4>(z,z,6) - l ♦ ^ (z,z)6 n - $ (0 *(z,z) + <$<fr^Cz, 2 ) ... (3.52) 

n*0 

If we proceed one step further by expanding p(a) in power series around s«z. 


then we have 


Art# 
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P<i>- l rt (,) (I| - P<») + ... 

rO 


(3.53) 


where 


r» d p — 

0 ■ «■" • 1 ■ 6x 


with 6<<1, and total derivative with reapect to z can be written at a sum of 

* 

partial derivatives , 


-^--^■+54= 

dr 3z 0 3z 


(3.54) 


We substitute eqs.(3.54), (3.53) and (3.52) into eq. (3.50). Using 
the fact that p(s) is slowly varying over the support of It, it can be taken 
in front of the integral as p(z). Thus eq. (3.50) takes the following form 

E <^ n \z,z)6 n -j l 6 m4 ‘ 1 p^ m \'i')fdBf dz" ~ y _1 Xy l <W z '', z ")6 n }rect d (z) 

n«0 lm-0 * * n«0 ) 

We also need to expand ^ n ^(z'’,z 1 ') around z"“z. Doing that we get 


4> <n) (z',z") * '!#’’,z) + 6 |4(z',z) (z'-z) + ... 

Using eq. (3.56) in (3.55) and by equating coefficients of $, we find 

<S * 0 * : fz ♦*°*( z »*> " 0 

( 0 ), 


therefore f {z,z)«4i v '(z) 


(3.55) 

(3.56) 

(3.57a) 

(3.57b) 


dz ♦ < ‘ 1) ( z » z ) + f^z.z) « p(z) | |ds|dz > y“ 1 (z)xy(z")A*'»*)jrect d <j5) 

' * /•% CO 


(3.58a) 


From eq. (3. 58a), we obtain by using (3.57b) 

z 


^ (1) (z,z)--^^z + p(z) ||dz'|ds|dz"V 1 (z")yy(z^)K®fe)rect d (z)+^ (1 U) 

0 ' (3.58b) 

Following the two-variable procedure, which requires that the term does 
not grow as fast as z, l.e., 
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lim 1 ^ CD ^ ^ „ 0 

z-h» z y ' ' ' 
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(3.59) 


In Appendix B, the condition given by eq. (3«5S) ie evaluated. The recult 
lead* to the following equation for ^^(z): 


p(*)M^°^(*)rect d (*) 


(3.60) 


and employing 


- k ■ ■' 0 — 6 ' 0—Q » 


(3.61a) 


6 eft-,-} 

0.0 . 0 0 . n.o 


(3.61b) 


If we use eq. (2.26) along with the fact that h.°h° + v^v++ In 


eq. (3.61b) , we have 


m 6 - &1± f(i°,i°) 


(3.62) 


Equation (3.60) can now be written as; 


d ,0 2 tt1 — - 0 . 

S* (=)-- — P(z) f<P (z) 


(3.63) 


♦ ( 0 ) (z) - 


l (i ° ,i° ) 0 


o -I yf.i?) 


(3.64) 
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Equation (3.63) represents two Independent equations, namely 


dA°<«) 


dz 


~r *<*> I .1“ )-a°(T) 


(3.65) 


I ( i+ 

dz 

The equations (3»&5) and (3.66) are two first order systems of 
differential equations for slowly varying coefficie s. 


(3.66) 


C. The solution of the slowly varying equations: 

Here we will solve the slowly varying differential equations in order 
to obtain the mean field in the air, the ground and the vegetated slab as shown 
in fig. 3. We assume the transverse mean field is excited by a plane wave 


- 0 " iK o z 

<'*' ± (z)> * y\a_ e 

qe{h,v} , 2 > 0 

(3.67) 

having polarization ^ and propagation constant 

V k o CO80 o* 


We can write the mean field in the three 

regions using eq. 

(3.32) 

and eq. (3.67). It is: 



w - 0 -itc_z , + 0 iic-jz 

<¥(z)> “ y*q e 0 + y*T «q e 0 

z*0 

(3.68) 

<H'(z)> ■ y*A(z)e iK 0 z + y^B(z)e* K 0 Z 

0>z>-d 

.(21.69) 

„ _ -i< z 

<¥(z)> «y *A e 
g -g 

z<-d 

(3.70) 


where T i« the dyadic reflection coefficient and A is constant vector. 

"3 ' 8 

The subzerc on quantities such as k. and k which indicate an angle of 

u Sq 

incidence 6 q will be suppressed through the rest of section C for nota- 
tion convenience, and tc e is given by 

B 0 

k b * k n< e ~ ein 2 0 ) 1/2 (3.71) 

g 0 0 R 0 

The y and y + in eqs. (3.68)— (3.70) are given by eq. (3.36) and A(z) , 
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B(z) are the slowly varying coefficients. Assuming constant P, from eqs. 
(3.66) and (3.67), we find 


A(z)«e " ‘Aq , 

where A. and B~ are constant vectors. 
-0 —0 

can be expressed as_ 



where 


i~-p f (ij.i?)* 

B(z)-e K " '*<, (3-72) 

The mean field eqs. (3.68)-(3.70) 

z>0 

0>z>-d (3.73) 

z<-d 



A(z) 

B(z) 



(3.74) 


Boundary Conditions: 

(i) At z-lT-'O, the tangential field has to be continuous, hence 

r<y,> - r<$ > (3.75) 

X S 

where <¥*>«y$ and <Y > m y<p represent the mean fields for z>0 
r 8 8 s 

and 0>z>-d, respectively. 

After substituting eq. (3.73) into (3.75), we multiply both 
sides by y 1 . Then using eq. (3.43), we obtain 



z*z“0 


(3.‘76) 


Using eqs. (3.74) in (3.76) we conclude that 

q° - A(0) - Aq (3.77) 

r«q° - B(0) » Bq 


(3.78) 


where 

and 

where ' 
ground 


where 
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(il) At *■- d, sT— "i, the tangential field has to be continuous, hence 


r<f > - r <$ > (3.79) 

* 8 

where <*? >- y<fr tt represents the mean field for *<-d. 

O ft o 

After substituting eq. (3.73) in (3.79), we multiply both 
* -1 

sides by y . Then using (3.43), we obtain 

* 8 (-d) - (y^ry^g (3.80) 

-1 

Substituting <p , 4 and (y Ty ) into eq. (3.80), we obtain 
8 8 8 

(3.81) 

(3.82) 

(3.83) 
(3.84) 

. and r are dyadic transmission and reflection coefficients at the 
respectively. 

Using eqs. (3.83), (3.84) in eqs. (3.81) and (3.82) we obtain 


, T _ -i(K— K )d 

A( -d) - ^H(y)ry g }-Ae 


1 + ; t _ i(K+K )d 

B(-d) * 2^-tCy )fy^A,e 8 


8 6 
A. 

tiy, m r\K/V.Q and y~ is given by eq. 3.36 for z<-d. 
Now we Introduce 
A - T-A(-d) 


B(-d) - T.A(-d) 


1(*C”K )d 

j ■ ^ e ® 

- -8 


T 

-8 



h°h° + T 


8 V 



(3.85) 


(3.86) 






where 
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(3.89) 


In eqs. (3.89) r h and r v ere Freenel reflection coefficients which 

5 ** o v 

are given by 


K ->*_ 


-& 


gh K + K 


gv 


e g K - K g 

e k + k 

8 8 


(3.90) 


We also obtain using eqs. (3.72), (3.77) and (3.78) 


k d Ik d 

r ■ e • r «e 

■8 «g 


(3.91) 


where 


. K I± + 2IP T (g'gy 


(3.92) 


For the final solution, we write the mean field in three regions as 



r o 


f _ 0 1 


q_ 

-itez , 
e + 


<n*,q)> - 



. ri(i 0 _xi“>a°_ 


,h(i+*i + )£ 8 *a? - 


Ikz 


, z>0 

(3.93) 


and 



T 


r + i 


A. 

-ijfz o 

i. 

as 

<'f(z,q)> - 

n(i°_ xi") 

m U , 

* e ‘3„ + 

_n(4*l + >_ 


and 


<$(z,q)> ■ 


■g 


lVA. 


lie d q -itc (z+d) 
e “ • e 8 


2 0 

• e • T »q , 0>z>-d 

-S-- * 

(3.94) 


, z<-d (3.95) 


where ig is i®. (z) for z<-d as given in eq. (3.37a). 

The first term in eq. (3.93) is the Incoming wave and the second term 


represents the reflected wave after travelling through the equivalent medium 
and having been reflected from underlying ground. The same is true inside 
the equivalent medium except that propagation constant now is k as defined by 
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eq. (3.92). The original wave after decaying through the equivalent medium 
is partially transmitted into the ground where it travels with a propagation 
constant k . These results, in the case of low frequency (Rayleigh regime), 
are identical to those obtained by Lang and Sldhu [1982]. 

Transverse Dyadic Green's Function in Three Media 
To compute the backscatterlng coefficients we have to evaluate the 
transverse spectral density. To evaluate it we will need the transverse 
Green's function in the equivalent medium which gives the scattered field at 
x in response to a unit dipole located at x ' . The Green's function will be 
calculated in the equivalent medium and thus, is an averaped quantity. We 
have defined the Green's function problem at the end of chapter II. Now, if 
we take the transverse Fourier transform of eq. (2.57) and proceed in the 
same way we did in the mean wave case, we obtain 


(f 

LG^.z, z') ■ I6(z-z')+ ids dz^ > p(s)x(k t »Z“S,z'" , -s)G(]^.,z'■*,z' , ) rectj(z) 


+ radiation condition 


(3.96) 


The definition of x 7 , I and £ are given by eqs. (3.16), (2.58) and (3.14), 
respectively. 

We write 

GQ^.z.z') - G(z.z') - y(z)G(z,z") (3.97) 

i 

Substituting eq. (3,97) into eq. (3.96) and using the Identical pro- 
cedure as we did in the mean wave case, we find 


£g(z,z') - y“ 1 (2r)6(z-z")+|jdsJdz"p(s)y“ 1 (z)t(z-s,z"-s)y(z-)G(z-,z')| 


• rect £ j(z) 


(3.98) 
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We would like to put eq, (3.98) Into a more familiar form by multi- 
plying from the right aide by Ty(0, and using y'^Ty-I^ we obtain 

ffit.z') - I^(z-i')+|jd.|<) t "c(.)y' 1 (i)t<«-..«"-*)y(«")GU" > r')J 

• rect d (z) (3.99) 

where 

A 

G(z,z') - 6(z,z')Ty(z') (3.100a) 


From eq. (3.100a) we have 

G(z,z') - G(z,z‘')y“ 1 (z') (3.100b) 

Using eq. (3.100b) in eq. (3. $7) we obtain that the transverse Fourier 
transform of Green's function is in the form of similarity transformation, 
i.e. 

G(z,z") - y(z)&(z,z')y~ 1 (z') (3.100c) 

We will solve eq. (3.99) approximately by utilizing the two variable method. 

We expand the Green's function in a power series and assume that the unit 
source is inside the slab, i.e. 0>z^>-d. We have 


G(z,z') - G(z,z;z‘*,z",6) - l G (n) (z,z f z'\z')6 n 

_ _ n«0 

where z ■ 6z and z'»6z" with 6«1. 

Proceeding in the same manner as in the mean wave, we obtain t 

% 

6(0>: - I^d(z-z') (3.101a) 

From eq. (3.101a) we have for zj*z' 


G^(z,z;z',z') - &*°^(z;z',z') 

+ G (0) U;*'.*O-F(*)||ds|dz > 'y" 3 ^jG (0) (z;z%z-)rect d (z) 


(3.101b) 


. a***- 
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Following the two variable procedure which requires that the G 
term does not grow as fast as z, i.e., 

£ i ft' 1 - o 

then we obtain 


^G.-^PfG 


( 1 ) 


(3.102) 


where f is given by eq. (3. 64) and 
G * G^(z,z"',z") 

A 

In eq. (3.102), 6 represents slowly varying terms for three media, and 


in each medium is noted by 


e + 

s 


z>0 

0>z>z" 

z'>z>~d 

z<-d 


(3.103a) 


with 


G f- 


■ 0 

0 * 

■A . 

V (z) 
“ 8 1 

A * (z)1 
s 2 

A 

fA (z) 
“g x 

mo (Z > 
g 2 

(7) 

L i 

(7) 
2 . 

’ G s- 

2s <«> 
» 1 

tr (z) 

“ S 2 J 

• V 

0 

m 

0 


(3.103b) 

I 

Now the boundary conditions to be used are that I'G(z,z'',6), has to be 
continuous at z«0 and z— d, In addition. Green's function must satisfy the 
jump condition at the source. 

The boundary condition at z-0 gives us 


A (0) - 0 , i-1,2 


(3.104) 




and 
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i i 

Th« boundary condition at z--d, z-~d gives us 
- - -1 T _ K< -*c>d 

V” • 4 [y ?v • v * «•* 

_ i J . l(K+»)d 

v d) - 4 [y r v • v 1-1 • 2 


(3.105) 


(3.106) 

(3.107) 


Equations (3.106) and (3.107) are the same as the ones we obtained in the mean 
wave case. Proceeding in the same way 


i g . - Si-£, ( - d > 


and 


b“ c-d) « r -a“ (-d) 

■s . «i *s . ' 

1 1 


i~ 1,2 


i-1,2 


(3.108) 


(3.109) 


Therefore following Appendix C, we obtain 

i(K-Od 


T„ « T - T ■ T e 


8 


41 *2 


“8 


(3.110) 


r. - r - r - t e 

1 2 “8 


21<d 


(3.111) 


where T and F are given in eqs. (3.86) and (3.89) respectively . 

“8 

Integrating eq. (3.101a) from z'>e to z'+e and letting e-+0 yields the 

t 

jump condition across the delta function. It ie: 


£<*'> - 


(3.112) 


From eq. (3.112) we have 


A + <z')-A~ (z') - J A 
*8 1 ' ' -s^ ' -A i 


1 * 1,2 


(3.113) 
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and 


B + (*') - B" (z') - 

" 8 i -i " B i 


i-1,2 


(3.114) 


where 


J A - 1 

" A 1 " 


V° 


(3.115) 


j B • 0 
B 1 


h -i + 

B 2 


(3.116) 


We obtained fromeqs. (3.102) and (3.103) 


+ _ + _ ifr? + 

».<*>"• ' io , I #,<*>■• - So . 


, 1-1,2 z"<z<0 

(3.117) 


. _ . 

A (z)-e -A 

S i *°i 

where l ±m T(i® ,i° + ) . 


,2itp -?♦- - 

- - 1 ir £ z - 

B (z)-e »B 

-s^ =0 i 


, 1-1,2 -d<z<z 


We assume Aq^ and are known, then the rest of the terms will be ex- 
pressed in terms of AJT and B+ . After that, using jump condition A~ and B* 

u i °i -0^ -Ot 

will be determined. We have: 


. _ . 

A (z)-e • A 

-s i N ' -0 1 


. _ Mfz . 

- r -e\ • 1-1 ■ 2 


(3.118) 


i^r d 

A ■ T;e K ’ . a " 

•gj -1 -Oj 


- t •' 

i 0 •- -J A 

u i A 1 


, i-1,2 


(3.119) 


and 


R + _ . tc 

— e 

-Oi 


-illl y i- 


• j „ + r . a „ 

-B^ -s -0 i 


, i-1,2 


(3.120) 


Where J A , J_ and T are defined by eqs. (3.115), (3.116) and (3.110), 
A i B 1 1 

respectively. 
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We also obtain for £ fi the following expression 


i* + d 


1< d 


r - e 
■8 


I. • 


1 - 1,2 


(3.121) 


where tc and £ are given by eqs. (3.92) and (3.89). 

© 

Substituting A's and B'a Into the eqs. (3.103) and using eq. (3.100c) 
we can write the Green's functions In three aedla, It will only be 
needed in the region 0>z>z". 


G s (kt ,z,z') - y(z)S 8 (e,z'*)y~ 1 (z‘') 


(3.122) 


where the slowly varying tern 6 Is given by 


2+ 

G s < z.zO- 


.2irp7+-— 

± ~Ti * + 

e *B 


i^£f + z . 
kC m r 
e -g 

0„ 


(3.123) 


and B n and B n are defined by eq. (3.120). 

U 1 2 

—1 2 

Substituting y(z), y (z") and G_ Into eq. (3.123), we obtain 




*11 


*21 


*12 


'22 


(3.124) 


whe re 


£u 


““o ( *S + * is" *' 

"2k I e +e ' Is '' 


(3.125) 


Sl 2 ' wj e 


i£ + <*-*'> „0 . i£ +I . . i- ~ * 


( iTxI ) + e 


0 


£ * e .(1>I )V (3.126) 


cot COt 
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‘o(„o T +, ‘fiV*') 0 ,+. ‘s +1 „ 1 s" 1 

_»/4 ~T V- + (jLf X | )"• *Ig*« 


.J^L 0 +. 1 

21 2k * 


-!!o 

22“ 2k 


with 


S O + 1 £ + <*“*') 0 + 0 + **”* n - ) 

Oj x i >• . (i;xJ + )+(lJxi + >* •£.•• • 02xl )| 


0>2>*“>-d 


G s ( x,jc') 


G s (k t .*.*-')e 





ORIGMt 
OF POOR QUAUTV 



(3.127) 


(3.128) 


CHAPTER IV 



ELECTROMAGNETIC EACKSCATTERING COEFFICIENTS FROM 
“ A LAYER OF VEGETATION 


In this chapter we will illustrate the application of the methods 
developed in the previous chapters by calculating the backscattering co- 
efficients from a vegetated slab of scatterers that are comparable in size 
to a wavelength. The physical configuration is shown in fig. 2. 

First, we will evaluate the transverse Fourier transform of the 
scattered field. The knowledge of the transverse mean field due to an in~ 
cldent plane wave and the transverse Green's function in the equivalent 
medium, will be used in conjunction with the spheral density to obtain the 
general form of backscattering coefficients. 

Then we will obtain explicit expressions of the electromagnetic back- 
scattering coefficients in terms of the bistatic scattering cross section of 
an individual scatterer. After determining the backscattering cross sections, 
their behavior will be studied by associating each term in the total back- 
scattering cross section with a physically understandable scattering process. 

i 

The Transverse Fourier Transform of the Scattered Field 
In chapter II we found, by employing the distorted Born approximation, 
that the correlation of the fluctuating component of the scattered field is: 

<V £ (*)ff(i)> “ |dsp(s)4' (x.s)^ (x,s> (A. 1) 

~ V e e 

where V is the volume of the slab of vegetation and 4* (x,s) is the scattered 

8 e 

field at x In the effective medium. It la given by 
“ 49 

?•'* 


lif 

* 
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Yg (*.£) - [ d *.'G (x , x' ) | dx' 't '-a ) <f (x" ) > 


From chapter II, we have 


*<*»*') ■ 


Su 


(4.2) 


(4.3) 


Bith £11 * i <*•*'> 


and t(x,x") 


(2tt ) ' 


dk dk' t(k,k")e i( - - “ 


(4.4) 


We have written in chapter III for the mean wave 

•A » A 

<4'(x'")> - <nk ,z") e , k -k.sine^ 

-t 0 ~ t o 0 


where 


<V(k. ,z")> - y(z'"*)4 , !C ,( z ^ > ) 
~*0 8 


(4.5) 


(4.6) 


The Green's function used in eq. (4.2) can be written as 

G(x,x") * G(x -x' t z,z') (4.7) 

— — t t 

The Green '8 function can be represented by its transverse Fourier transform: 


GC^ 


-x'iZ, z') - — fdk^. G(k 4 .,z,z')e 
“* (2 tt) Z i 


%t* ( *t ’5t ) 


-t '^t’ 


(4.'8) 


In order to obtain the transverse tourier transform ox eq. (4.2), we 
substitute eqs. (4.5) and (4.8) into eq. (4.2), and also use 

i(kXx; * k^xD 

j |dk;dk;"t(k;,z';k;-,z--)< 

(2 W 


t(*",x~) - — [dkrdkr"t(k;,z';k;%z^)< 


(4.9) 


In eq. (4.9) t(k',z';k'',z") is different from £<k',k"). We have used the 
same symbol for convenience of notation. 
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Now we have 


i (k^.z.s) - (2n) 2 fdz"d*^G(l^,a,Ot(k <r ,*"-iil^ ,*"-•)<$(]*,. ,*"): 

8 * “* J -t -t -t 0 




(4.10) 


From eq. (3.122) t we have 

GQc^z.z') « y(jc t ,z)G(k t ,¥ f ¥')y“ 1 (k t ,z'‘) (4.11) 

-1 * 

where we have explicitly shown the dependence on 1^. in y, y and G. We now 

use eqs.(4.6) and (4.il) In eq.(4.10). The Greens function GQc^z.z") and the 

mean field 4> (k. ,z ") are replaced by G(k ,z,s) and <fc (k ,s). This is a good 
s ”t Q c s c 0 

m 

approximation to first order in 6 due to the behavior of t about z' * s and 
z" ■ s. The result is: 


c )-s 

* B (k t * z »s.) “ f g (k t ,z»s)e 0 

e e 


(4.12a) 


where 


4*3 (k^z^s) - (2Tr) y(k t ,z)G(k t ,z,s)p(k t ,k t ^, 8)^(1^^, s) 


with s • 6s and 


(4.12b) 


PCJ^t*^* 8 ) * | dz "d z 'V 1 (k t ,z")t(Jk t ,z"-B;k t ^,z > '-s)y(k t ^,z^) (4.13) 


The Formulation of the General Form of the 
Backscattering Coefficients 

The backscattering coefficients are directly related to the transverse 
Fourier transform of eq. (4.1) with respect to and x* evaluated at the 
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upper interface (z-0). We atart by taking the transverse Fourier transform 
with respect to 3^ and & t of eq. (4.1). We have 


<? f (k t »*)¥j<£ t »*) > ■ pfdsfa (it****) 

i ( e e 


(4.14) 


where p is the constant particle density. 

Now, by putting eq. (4.12) into eq. (4.14) and by integrating over 

•k 

and by setting z*z*0, we have 


<v £ (^,0)? + (k t ,0)>-s(i^,q|ii)fi(iCt-it) 

where q 

S(k ,q|u>) - p] dst (k t ,0,B)f^ (kj.,0,8) 
id e e 


q e {h,v} (4.15) 


(4.16) 


Here SOc^qjos) is the transverse spectral density of ^(x) at the interface. 
The symbol q represents the polarization of the incident wave and u the orien- 
tation of a typical scatterer. The spectral density S is a 2x2 matrix of 
dyadics . We write 


S(k ,q|u>) 


11 

«u 

21 

*22 


(4.17) 


where and are the dyadic spectral densities of electric and magnetic 
fields respectively. In our development, we will concern ourselves only with S^. 

It is shown by Lang [1981] that (Lang has used S for S^^) is directly 

0 1 
related to the backscattering coefficients o by 

pq 


2 2 

_ k^cos 6„ „ 

O 0 _.^__ 0 £ 0.3 q|u) . £ 0 


pq 


p,q e {h,v} 


(4.18) 


4tt' 


where 0 Q is the angle of incidence of the plane wave with respect to the slab 
normal, and pP,q^e{h®,v®} ate unit vectors indicating the polarization of the 
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To simplify the remaining computation for S^, we note from eq. (4.18) 


that S(k t ,q|u) will only be required for kg ■ -kg 
verse transition operator t in eq. (4.13) cau be expressed 


-kg a in 6 qX . The trans- 



_ 1 _ 

2tt 



,k , K - )c ik'(r-.)-iK"(z 
”*0 




(4. 19 > 

Using eq. (4.19), y"* and y from chapter III, we compute eq. (4.13); 
then we substitute this value of p , y, y , G and $ g fro® chapter III in 
eq. (4.16). We find that the spectral density is: 

4 ? 4 

^(-kj ,q| U ) - l (W t )(U*)d. (4.20) 


where 


-% 8 e “V 0 

W x - e • fj-e . S_ 

. + . + 

iK s k s n 

A aif) niQ 0 

-2 ■ £jj* ' ifi* * £|i- 

(4.21) 

“3 ' • lift' 4 s° 

-itc s ijc+s n 

w 4 - e * • l£_ 


with 



h m 

hi * * < 'VVV' 0> 

(4.22) 

\ 

hu * f<-\*- | b?*t 0 '- K b> 

“iv ■ 

(4.23) 


In obtaining eqs.(4.20) through (4.23), we have used the relationship between 

i + 

f and t given in eq.(2.26). In addition, k q and T are k" (as given in eq. 

8 0 

(3.92)) and T (as given in eq.(3.91)) respectively with k. replaced by k . 

* l£'<i - i£ + d . 0 0 0 0 0 

Finally 1.^-1.^- » . £g .. and £g -r gh A + T gy v. v + . 


■A 5 
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Substituting eq. (4.20) into (4.18) we have the general form of back,- 
scattering cross section as 



(4.24) 


Now substituting eq. (4.21) into eq. (4.24), we see that there will 
be sixteen exponential integral terms to evaluate. We examine these terms 
and find that phase cancellation gives 6^ order terms and phase accumula- 
tion gives 6^ order terms. 

Omitting the 6^ order terms, which are much smaller than 6^, eq. 


(4.24) can be written as 


0 0 . 0,0 
a * a . + a +o . 
pq pqd pqr pqdr 


(4.25) 


where 


°pqd * 4 ”° 




- - _* . 

-Iks -itc s „ r\ Iks. Ik. s n 
^0 r “0 o, 0, -n -0 0 

* q )p <e * f «e • q 


*(e f je 3 _)£_<e 


£)]ds} (4.26) 


- [f . . 1<„S 1 k 8 a -IK* 8 -1< S A 

j - 4Trp!][p *( r • e f * T e r.g ) p^.( I 1 * *e • f^*e r*q )]ds 
pqr v -s 0 -II “ £ ~ “ s 0 “ XI 

**d 


+* +* 
-Ik* s -i< s 


(4.27) 


and 


0 / (fir 0, -%% 1 £0 8 _ 0, 0. 1 S0 8 * _i -0 0. . 

r pqdr " 47rp jj j [ £- (e * -iv e * I^iJE < e * liv e c 

— d ® 


, . o , ; i 5o% “^o 8 o. o , :* “ 1 -o 8 * i 2o s o.. 

+ [£-•( I 8() - e • Irif • 3->£- ,( Is* e * Inf • S_>1 
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^ r o, -i £b B - o. o . -i £o *.* 1 -o • o* , 

+ [£_•(« • f IV *« • r.jj 2 . • < £^j • • • j_)l 




+* 


a / i£ n 8 ~i<-8 n A i£rt • * -iic_ a . n ) 

+ "°q_)fj(e * | IV Pe • f-aJ]j d S| 


p ,q e{h ,v } 


(4.28) 


Here Opjjj, a p q r and °pqdr represent direct, reflected, and direct- 
reflected backscattering contributions as shown in figs. (4a), (4b) and (4c- 
4d), respectively. 


The Calculation of Backscattering Coefficients 

In this section, we simplify the expressions for the backscattering 
coefficients by assuming the scatterers have no average depolarization at the 
level of the mean, i.e., 


£ * £ 
hv + v + h 


(4.29) 


It was shown by Lang, 1981 that dipole discs having a uniform distribution in 
the 4> coordinate (azimuth symmetry) obey eq. (4.29). It is conjectured that 
eq. (4.29) holds for all scatterers having azimuthal symmetry. It should; be 
noted the eq. (4.29) also Implies that the average or bulk medium is uniaxial. 

Using eq. (4.29), the scattering amplitude can be written as 


~z ± ~z ± tO.O .-v± o o 

f ■ f. . h h + f vj , 
■ hh w — ±— ± 


(4.30) 


where 


+ 0 .,.0 .0, 0 
r ,t ■ *±-l<Af±±> * *± 


p e{h,v} 


(4.31) 
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Next a simplified expression for jcJ can be obtained by employing eq. (4.30) in 
eq. (3.92). It is 

*\ A- A #fc. 

(4.32) 

where &>* 0 1 


± ± v0. 0 . ± 0 0 

£o ,K k W +K v^±^± 


± , 2ffp — ± 

c ■ K« + ~A~ f 


0 K 0 pp 


“o * k o* ln9 o 


(*.33) 


Since Kq has no off diagonal terns, the exponential terns appearing In eqs. 
(4„26)-(4.28) can be written as 

* 4 * lK l z .0.0 . iK > 0 0 

e - e h h + e v^v + 

A. Calculation of like polarization cross sections: 

From eq. (4.24) we have 


(4.34) 


,° 

PP 


- u 4 

4tt P |e° J ± | i (W i )<Wj)|f°ds 


p e{h ,v } 


(4.35) 


where 


PP 


0 , 0.0 
a , + a + a , 
ppd ppr ppdr 


(4.36) 


We proceed by dotting both sides of eqs. (4.26), (4.27) and (4.28) 
with q^, using eq. (4.34) and taking the average over all particle orienta- 
tions gives 


°pp d - *«pIe3b“I 


-4Iidk d 

0 £ ._0i 2 1 - e P 


4Imtc 


p e{h,v) 


(4.37) 


'° " Attp |e°-xxEXI + 


.0.2 e 


4Imx + d 
P 


ppr 


4XmK 


- l p sP l 4 


(4.38) 


The Opp^j. has three components: 
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°ppdr 


0 

a ppdr x 




(4.39) 


where 


ppdtj" 


2In(K + -tc"’)d 
0,2 e p p - 1 


2Ib(/-k“) 
P P 


l r 


•P’ 


(4.40) 


a ppdr« “ 


0,2 e 


2Im(K + -K“)d 
P P 


- 1 


2Im(K + -K~) 
P P 


r spl 


(4.41) 


and 


PPdr, 


2Im(K + -K )d 

8irp Re [(PgjftP -MP^pJ) +~~ I r «r.| 2 




2 In. (/-<’) 
P P 


sp 1 


(4.42) 


with 


sp 


i(<I + K )d 
r e p P , 

8P * 


p e{h,v) 


(4.43) 


We note that If k* - k~ then eq. (4.40)~(4.42) change® as follows : 


2Im(tT - Od 
P P 


- 1 


2Iq(k - K ) 
P P 


(4.44) 


The Fresnel reflection coefficients appearing In eq.(4.43) are given by 


r - K °- l K & 

*h K fl +K g 


e < n -k 

* m _S _2 & 

gv C gK() +K g 


(«.«) 
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and 


*6 " k o CO80 o 


K g " W ln e o ) 


( 4 . 46 ) 

( 4 . 47 ) 


The Oppjj represents the direct backscatterlng contribution as shown in 
fig. (4a). The incoming wave propagates into the vegetation and is scattered 
directly back to the observer. When summed over all scatterers, eq. (4.37) 
results. The backscatterlng coefficient °p pr given by eq. (4.38) represents 
the sum of all waves which are first reflected from the ground, then scatter- 
ed, and finally reflected again by the ground towards the observer as shown in 

i 1 4 

fig. (4b). The reflection coefficient T appears in the equation since 

s p 

the wave In this case has been reflected by the ground twice. 

The third term 0 p pdr as given by eq. (4.39), results from two differ- 
ent but similar mechanisms as shown in fig. (4c) and fig. (4d). In one case, 
the wave is scattered and then reflected toward the observer; whereas in the 
second case, the wave is first reflected from ground interface and then 
scattered toward the observer. The third term, eq. (4.42), occurs due to a 

combination of both cases. These terms interfere coherently. The 
2 

|r | factor represents single reflection from the ground. 

©* . , 

0 0 2 0 Q 2 t 

The coefficients |p_.f|p_| and |p+fllEt+' represent the backscatter- 
. n 0 

ing cross sections in j>” andj> + directions, respectively. The coefficients 
Ije 5L I 2 and ljEL+ *!iv*E^J 2 re P re sent the bistatic cross sections in 

different directions. Note that the bar over the expression represents the 
average over angular variables. 
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B. Calculation of cross polarization cross sections: 

In this case the incoming and reflected waves have different polari- 
zations. Thus the equation for Opq dr cannot be reduced to a simple term aa 
in the previous two cases. But otherwise, following a similar procedure, we 
obtain 


0 0 0 . 0.0 
o' ■ a ■ a . + a + o . » p rq 

qp pq qp<» qp* qpdr 


(4.48) 


where 


V, ■ 


-2Im(»c +k )d 

0* 0,2 1 - e q P 


2Im<(c”+K ) 

q P 


(4.49) 


a “ 4 ttp Prf T * T q . 
qpr 


0i 2 e 


2Im(tc + +K + )d 

q p' 


2Im<K++Kp 


— |r | 2 |r | 2 

1 sq 1 ' sp' 


The a pq dr has three componentsi 


(4.50) 


a° «a° + o° f o° „ 

qpdr qpdr x qpdr 2 qpdr 3 


(4.51) 


where 


-2Im(K~-tc*)d 

a°_ - 4 tt P |,W|* 1_L_. |r I' 


' 4,p itaxrfi 


2Iiii(k"-k*) 

q p 


•P' 


(4.52) 


2Im((c + -(c")d 


J • 47 TP p .f • q 

qpdr, * p 'i-siVif 


W>i 2 q . p 


«■<<-? 


i f .,i 


(4.53) 
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8«p R« ^E^iiiS-XPiliv^^ 


* 1 - e q 

p r — - 

lq ip 

H * i(-K +K *1 

» n n 




K f K ) 

q p q p 7 


]i 


(4.54) 


+ - - + + — 

Note that If K p - K p or K q - K p + - K p =0 then one obtains the correct 

form of eqs. (4.52)-(4.54) by taking the limit as those equalities are 
approached. 


It is interesting to note that only the' imaginary part of the effective 
permittivity enters the expression for a pq . This is a result of the fact that 
o p q represents an incoherent sum of the scattered power from individual par- 
ticles embedded in the effective medium. 

If the particles 1 sizes are small in comparison to the wavelength, 
these results reduce to the results obtained by Lang and Sidhu [1983]. 
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CHAPTER V 


NUMERICAL RESULTS AND DISCUSSION 


Our main Interest In this chapter is to obtain the numerical results 
for the theory we have developed In the previous chapters. For this purpose, 
we use our method to model a forest canopy by a collection of lossy dielectric 
discs, which are assumed to have radius a, thickness T and relative dielectric 
constant Our formulation is applicable for arbitrarily shaped scatterers 

in the Rayleigh, resonant and geometric optics regimes; however, because of 
the availability of Rayleigh and geometric optic scattering amplitude algorithms 
we will limit our calculation to those regimes. 


Using the parameters encountered in active remote sensing of vegetation 
layers, we will present the numerical results for the skin depth and for the 
backscattering cross sections in the Rayleigh and geometric optics regimes. 
Before proceeding with this calculation, we will discuss the relative dielectric 


constant of the leaves and the scattering albedo of scatterers. 

The Relative Dielectric Constant of the Leaf 
Our calculation of the relative dielectric constant follows the model of 
de Loor [1968], and Fung and Ulaby [1978], We conclude on the basis of mea- 
surements conducted by Broadhurst [1970] that this model is not suitable jfor 


* — j t.-.i — *i mi- j..* 

i. X C^UCULXCO L/CXUW X UU£ UUC UW ItO J. OiiUl C LU 


_ r s j /-i — 

u auCvuia tut cuuuuclxvc ^uuuixt/ 


losses within the leaves; however, we will wse this model below 1 GHz in order 
to recover results obtained by Lang and Sldhu [1983]. The dielectric formula is 
given for the real and the imaginary parts of the relative dielectric permittivity 


below. 


e -5.5 
eZ - 5.5 + ; 


l+(fx)' 


(5.1) 
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Frequency in GHz 


Fig. 5. Real and Imaginary part of dielectric constant vs frequency 




5.5) 


fr 

l+(fx) 2 


(5.2) 


with e “ Z* + The letter T is the relaxation time of water which de- 

r r r 

pends on the temperature and f is the frequency. At 20°C, fi is approximately 
equal to 1.85/A, where A is the free space wavelength in centimeters. The e 

m 

appearing in eqs. (5.1) and (5.2) is the relative macroscopic static permittiv- 
ity of a leaf and given by 

' 5 + 51 - 56 \ 

where V is the volume filling factor which lies between 0.1 and 0.6. For 
m 

illustrative purposes, we plot a real and imaginary part of as a function 
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of the frequency as shovm in fig. 5 with V * 0.3. 

We observe that e' - '* has a maximum around 15 GHz and e** is constant un- 
til 1 GHz, then it decreases monotonlcally. We also have observed that chang- 
ing parameters did not effect the overall shape of the curves. 


i 



The Albedo of a Single Leaf 

In this section we will calculate the albedo of a single dielectric 
disc. This is necessary due to the fact that our theory is limited to low 
albedo particles where absorptive loss is then the dominant mechanism. First, 
the dyadic scattering amplitude for a lossy dielectric scatterer is required. 

In the Rayleigh regime, [see Lang 1981], the dyadic scattering amplitude f is 
related to the polarizability, a, by 



p,qe{h,v} (5.4) 


Ishiraaru [1978] then gives expressions for the polarizability of a lossy 
dielectric disc. 

The dyadic scattering amplitude for a lossy dielectric disc In the 
geometric optics regimes has been calculated by LeVlne, et al. [1982]. An 
approximate expression for the dyadic scattering amplitude is obtained by 
considering a plane wave incident on an infinite dielectric slab. The inter- 
nal fields In the slab are calculated ex ctly. Then the equivalent sources 
generated by these fields in the region of the slab corresponding to the disc 
are u&ed to calculate the scattered field. The approximation requires that 
the leaf be many wavelengths in diameter and have a thickness, T, small com- 
pared to the diameter. 

From LeVine, et al.,we have the following expression for the dyadic 
scattering amplitude in the geometric optics regime: 
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where A*e r -1., and the Fourier trensfora of the cross sectional shape of 
the disc is given by 

f iv or 

S^) - jS(xpe d*£ (5.6a) 


Witn j 1 s x: e disk face 

S(*p ■ I 

( 0 , x* i disk face 

where the prime quantities are coordinates in the face of disk and ■ 

Here and 0 t are the projections of !_ and () of fig. 6 onto the 

face of the disk. In eq, (5.5). F is related to an integral of the induced 
charge in the infinite dielectric slab. It is given by LeVine, et al, (1982], 
For a circular disk of radius a eq. (5.6a) becomes 


S(v t ) * ~~ J 1 (v t a) (5.6b) 

where J.(c) is the Bessel function of first order and v -IvJ. 

X t ' — t ' 

From Ishimaru [1978], we have the following expressions for total. 


scattering, and absorption 

cross sections. 



• (s. 

^Im[f(i° s i 0 ;q )£ 0 ] 

, qe{h^ ,v } 

(5.7) 

„« . 1 

|f(g°,i°;q)| 2 dn 


(5.8) 


Air 


and 

a a q> " (5.9) > 

V 

where f_(0®,i®,q) is the vector scattering amplitude due to an incident plane 


wave with polarization q. These vector scattering amplitudes are related to 
the dyadic scattering amplitude by 
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f(0°,i°) - f(0°,i°jh)h 0 + £(Q°,i° }v)v 


(5.10) 


In addition, in eqa. (5. 7) — (5 ^ 9) , dfl la the differential aolid angle, dV is the 
differential volume element, and ECx") Is the electric field inside the par- 

,<q> 


tide. Using the total scattering crosa aection, cr iJ and the acatterlng 

Vi 

rCq) 


croaa aection, o^^we can define the acatterlng albedo by 


w< q) - 

*0 


r<q> 


(5.11) 


The total cross section represents the total power loss from the 

incident wave due to the scattering and absorption of the wave by the scatter. 
One can show that [Born and Wolf, 1964], 



Fig. 6. The geometric configurer 


£ leaf 



t * 
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qe{h,v) 


(5.12) 


t a s 

The geometric configuration of a leaf for an albedo calculation Is shown In 
fig. 6. 

We evaluate eqs. (5. 8 ) and (5.9) for the horizontal and vertical 
polarization cases. _ The results at low frequency (Rayleigh regime) are: 

(5.13) 


and 


where 


°i h) - vr v 


sln^S, 


, v f o Bin d_\ 

■ vr v ^°’ 8 0 + TTff) 


r (h) _ V i a 1 2 


s 


Sr W ! 


,<v> . l£!i± 2 L s \ + 

s 6n y 0 


2 2 

V - | ira^T 


sln2e o j 


(5.14) 


(5.15) 


(5.16) 


(5.17) 


e - 1 
r 


(5.18) 


In fig. 7 we plot eqs. (5.13)-(5.16) as a function of the frequency 
at zero incident angle with 3=7 cm, T=0.3 mm and V *0.1, Typical dimensions 
of a leaf having radii of one to several centimeters and thicknesses of 
tenths of a millimeter have been used. i 

At zero incident angle, using the above eqs. we conclude that 0 W =cr (v). 

a a 

cross sections Increase with the frequency. For large values of frequency a 


s - a s 


, in/ V v / _ 

and «cr a =a fl as seen In fig. 7. The curves show that absorption and scattering 


is greater than a . 

a 

Following the LeVine, et al [1982] work, we evaluate o and a for the 

a s 

horizontal and vertical polarization cases at zero incident angle in the geomet- 
trie optic regime. 
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The absorption and scattering cross sections are given by 

alnh(K.T) 


and 


where 


V - * |< 1 *. = 


«,) «"> 2 °Sq 

T S " *Q COB 6 0 


W^ q ^(e,(j>) - |slnc(0 T )e^ + + sinc(e“)e^J 


[W (q) (0 r> ^ r ) + w (q) (e t ,* t )3 


H\ + 0 


- 0 ,2 


with Sq - S(0)“Tra^, C - l^ c o T ^ 2 011(1 


0* - [K ± k Q 0 x ']7/2 , 0; - n° • 0° 


(5.19) 


(5.20) 


(5.21) 


(5.22) 


Quantities not explicitly defined in eqs. (5. 19) -(5. 22) are defined 

4 * 

in the LeVine, et al, reference. We note that the e~ are forward and back- 

q 

ward going wave amplitudes in the dielectric slab for polarization q. The 
propagation constant in the z direction inside the slab is ic ■ + itc^. 

Using the same parameters we used for the Rayleigh regime, i.e. a=7 cm, 
T» 0.3 mm, and 0.1, we plot CT q and as a function of the frequency at 
zero incident angle (fig. 8). As we see from this figure, at a low fre- 
quency of around 1 GHz both of the curves are linear; however as the frequency 


y l 

J 


I 


Frequency in GHz 


Fig. 8. Scattering and absorbtion cross sections vs frequency in geometric 
optics regime 
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m ■ 7 cm 
T * 0.3 ton 
V " °*1 

n 




f ■■■ >■ ♦ 

.1 .2 .4 .6 1.0 2.0 

t 

Frequency in GHz 


9. Albedo vs frequency for Rayleigh regime 
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increases the a curve becomes greater than the a curve. 

s a 

Now, substituting 0 and a into the eq, (5.11) a low albedo check is 

8 t 

made for both regimes. We plot the albedo as a function of the frequency in 
the Rayleigh regime and geometric optic regime, as shown in figs. (9) and 
(10) with the same parameters (a"7 cm, T«0.3 mm, V *0.1 and 0*0°). In 
these figures it was observed that albedo Increases almost linearly in the 
Rayleigh regime and monotonically in the geometric optic regime. The ex- 
pression for the geometric-optic albedo obtained by using eqs. (5.19) and 
(5.20) in eq. (5.11) can be substantially simplified when the phase variation 
of the internal field across the slab thickness is small (k T«l). The 
simplified expression for the albedo is 

1 


W 


o 1 + o /a 
a s 


(5.23) 


where 


a cos0_ 

a r 0 




(5.2A) 


From eq. (5.24), we conclude that the albedo is independent of the 
leaf radius, but very sensitive to the thickness, incident angle, and water 
content of the leaf. 


O 


The Skin Depth 

A knowledge of the skin depth gives a better understanding of the 
mean wave’s behavior inside the vegetation layer. The skin depth is defined 
by the depth of penetration at which the wave’s amplitude decreases to e~^ 
of its initial value, we write 


Skin Depth 


Im k 


pe{h,v } 


(5.25) 
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where from eq. (4.33) , ve heve 

K p - V 0,e o + 5^5^ 7 pp ( - } pc{h,v) <5l26) 

We heve used the feet the I . -7. -0 end k*-k”»k . These result* esn be obtelned 

vh hv p p p 

from (Leng» 1981) for Rayleigh discs end from Appendix D for geometrlc-optlc discs. 

We can see from eq, (5.25) that If the vegetation depth Is large com- 
pared to the skin depth, then the effect of the ground reflections are neg- 
ligible, whereas, when the skin depth Is large compared to layer thickness, 
the characteristics of the ground become Important In interpreting radar 
backscatter information. 

An examination of eq. (5.26) shows that the skin depth for both h 
and v polarizations contains a factor of cos0q in the numerator. This 
forces the skin depth to zero as the angle of incidence approaches grazing. 

Due to the increasing angle of Incidence, the effective depth of the propa- 
gation keeps decreasing monotonically as seen in the following figures. In 
fig. (11) and (12), the skin depth is plotted as a function of the incidence 
angle for both polarizations in Rayleigh and geometric optic regimes, re- 
spectively. In figs. (13) and (14) the skin depth is plotted as a function 
of the frequency for two different regimes. 

Now, let us examine figs. (11) and (12) in more detail. In the case 

of horizontal polarization, the electric field is approximately parallel to 

1 

the leaves at all angles of Incidence and thus the skin depth has onlycosQ^ 
dependence. The same thing is true for vertical polarization, but only when 
the angles of Incidence are small. Thus, both have the same skin depth at 
small angles. As 0^ becomes larger, the electric field tends to become per- 
pendicular to the leaves and each leaf absorbs less energy. This explains 
the increasing skin depth as a function of the incidence angle for vertical 
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polarization. Eventually, the skin depth decreases due to the cosd^ In 
the numerator in the vertical case. 

The curves of Rayleigh and geometric optics regimes have identical 
shape as seen from figs. (11), (12), (13) and (14). The only real differences 
are that the curves for Rayleigh regime decrease more rapidly with increasing 
frequency and the skin depth is much higher than values obtained for the geo- 
metrical optic regime. These results were expected due to the fact that the 
physical size of the scatterers become an Important factor in the high fre- 
quency regime. 

Both figs. (11) and (12) have the same parameter: e-7 cm, T-.5 mm and 

P"500/m . In figs. (13) and (14), we have choosen a“7 cm, A value ofT«.5mm was 

picked for fig. 13 to recover Lang and Sidhu’s result, but in fig. 14 we have 

used T*.l mm in order to obtain a low albedo. The Rayleigh results have used 

scatterers with polar angle 0 distributed uniformly, between 0 and A0||"3O° (see 

Lang, 1981) and the geometric optic results have used scatterers with a fixed 
0-30°. 

The Curves of the Backseat tering Coefficients 
We developed a computer program that provides numerical solutions for 
the backscattering cross sections which we have obtained in chapter IV. In 
the Rayleigh region the results are identical to those obtained by Lang and 

Sidhu[1983] . For purposes of completeness, we present some of their results in 
through figs. (15), (16), (17) and (18). These figures have common para- 
meters: a ,s 7.5cm, T“0.5mm, p*500/m J ’, f«400MHz, £ r *30. 8+10.62 and £ *12+13. 

© 

In fig. (.15), and ^hr are shown as a function of 6g* In 

this case, the term is greater than for all angles of incidence. 

This means that the energy backscattered from the ground is a significant 
part of the total backscattered energy. Since the skin depth is large 




Angle of Incidence (deg) 


Fig. 16. Backscatterlng coefficients vs angle of Incidence for Rayleigh 
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compared to the layer thickness for nil angles of Incidence , we are re- 
ceiving scattering back from the same number of scatterers for each angle of 
incidence, resulting in a flat behavior for The reflected terms 

and slopes up with an increasing angle of incidence with the increasing 

reflection constant. 

* 

Figure (16) shows the total backs cat terlng coefficients plotted for 


w 


a perpendicular distribution of leaves. At small angles of incidence a 

0 0 
and (Jjjk are equal, and as the angle of Incidence increases, a decreases 

because the vertically polarized wave becomes parallel to the leaves and thus 

0 

gets absorbed. On the other hand, has a flat response because the chang- 
ing angle of incidence does not affect the polarization orientation with 
respect to the leaves. 

Figure (17) shows the total backscattering coefficients plotted for 

a parallel distribution of leaves. In this case as the angle of incidence 

is increasing, the horizontal polarization is not affected since the leaves 

are parallel to the polarization at all angles. The backscattering cross 

section is higher as compared to the perpendicular distribution at small 

angles of incidence. Comparing o®. in figs. (17) and (16), we see that the 

vn 

vertically inclined leaves give rise to more depolarization as compared to 

the horizontally inclined leaves. 

To see the effect of moisture in the underlying ground, we have i 

selected two different dielectric constants. For dry ground, we have used 

e ■ 12+13 , while for very wet ground we have used the dielectric constant 
8 

of water obtained from Debye's formulation [1929] » i.e.. 
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with 


and 



(5.27) 


(5.28) 


where ^ is the wavelength in centimeters. From eqs. (5*27) and (5.28) we 
have Eg»80+il.5 at f*400 MHz. Figure 18 shows plotted for these two 
different ground conditions. From this figure we see that there is a large 
difference in the backseat tering coefficients for different ground moistures 
because the ground reflected terms are dominant. When the layer thickness 
is increased, the ground reflected terms become small. Hence their effect 
is masked and only the direct component, a p p( j» pe(h,v) is important. 


Next, we present the backscattering curves for the geometric optics 

regime. All of them have the same parameters: f*5 GHz, a* 7 cm, T»0.1 mm, 

and e “9.75+11.31. An example of the results is shown in fig. (19) for h 
r 

polarization. In this computation the discs all have an inclination angle 
of 30° with respect to the slab normal (0=30°), however they are uniformily 
distributed in the azimuth coordinate. To aid in the interpretation of the 
results, the direct, reflected, and direct-reflected components have also 

i 

been plotted. An examination of fig. (19) shows that has two major 
peaks: one at an Incidence angle of 0° and the other at 30°. The peak 

at 30° comes from the direct and reflected components and is due to specular 
reflection from the front and back faces of correctly oriented leaves. The 
larger term at 0* is caused by the direct-reflected components of the 
backscattering coefficient and is caused by the transmitted energy through 
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the leaves that are reflected from the ground. The width of each peak Is 
related to the radiation pattern or the beam-width of an Individual leaf. 

The leaves considered have been chosen to be thin. This was necessary to 
keep their albedo small so that the distorted Born results could be applied. 
It Is anticipated that by Introducing a distribution of leaves, which have 
an inclination angles that depart somewhat from 30*, a considerable smoothing 
of the curve will occur. The parameters which we used In fig. (19) are 

p«500/m 3 , d-1 m, 0-30*, and e -73.5+121.1. 

6 

We again obtain z from Debye’s formula to see the effect of moisture 
© 

In the underlying ground. We use the parameters of fig. (19) In order to 
plot fig. C20) with a new ground dielectric constant e -13+12. This value 

O 

corresponds to volumetric water content (cm^/cm^) of 0.3, which is relatively 
dry ground, Wang and Schmugge [1980]. 

From figs. (19) and (20), we can see that, except for small incident 
angles, there is almost no difference in the backscattering coefficients for 
different ground moistures. This result differs from the result obtain in 
the Rayleigh regime as is shown in fig. (18). The difference is due to the 
physical optic scatterers which we are employing. They have a gain in the 
direction of the Incident wave. It is this gain that was essentially com- 
pensating for the reflection loss at the interface. In fig. (20) we also 
plot and for illustrative purposes. i 

In fig. (21) we use the same parameters as in fig. (20) except p and 
d are increasing to values of p-2000/m 3 and d-10 m. We plot c^, o^ v and 
0® in this figure. A comparison of the two figures shows that they are the 
same except that fig. (21), no longer has the peak at 0- 0°. This is because 
increasing P and d made the skin depth substantially smaller. Thus the in- 
cident mean wave never reached the ground and there was no dr component in 
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Che backseat ter. The layer essentially looks like a half space. 

In fig. (22)., we illustrate Che effect of 6 on Che backscaccerlng 

cross secdon In Che geometric optics caee. Since 6*90* , this figure 

corresponds to the leaves perpendicular to the ground. We also plot o^ and 

o„t, on the same figure. In this case, we have cfi Jkfi a 0 and o® KsP . . 
vn • ppd ppr pp ppdr 

We see that has flat response due to the fact that the leaves act like 
a corner reflector. 

At zero angle of Incidence the value of and are equal, but as 
the angle of Incidence increases, decreases because the vertically po- 
larized wave gets absorbed. It was observed that figs. (16) end (22) follow 
almost the same trend, except $he depolarization is smaller in the geometric 
optics case. 

So far, we have shown some backscattering coefficient curves. These 
curves are not smooth as obtained in the Rayleigh regime. This is due to 
the radiation pattern of the scatterers. In the Rayleigh regime, the dipole 
radiation pattern has a large beamwidth which in turn gives a smooth response. 
But in the high frequency domain, the radiation pattern has a smaller beam- 
width compared to the Rayleigh regime, and the scatterers behave as a strong 
radiator. This in turn gives us the backscattering coefficients as we have 
shown throughout figs. (19)-(22). This is the reason why we observe a peak 
at 30* in the geometric optic regime. 

We have analyzed numerically discrete scattering model for vegetation. 
We have calculated the skin depth and backscattering coefficients explicitly 
in terms of the scattering amplitude. The scattering albedo was calculated 
to make sure that it is small. It was found that the leaves had a much 
higher albedo than was previously anticipated at microwave frequency above 


5 GHz. 
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It was found that for vegetation layers less than one skin depth 

thick and for moist ground that a° . was the dominant term, It was also 

pqdr 

observed that thin layers had ba eke catte ring coefficients that were almost 
flat functions of angle of incidence in the Rayleigh regime. In the 

geometric optics regime this was not the case because the scatterers which 
we were employing had gain in the direction of the incident wave. It was 
this gain that was essentially compensating for the reflection loss at the 
interface. 


c 
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APPENDIX A 

RELATIONSHIP BETWEEN THE TRANSITION OPERATOR 
- AND THE SCATTERING AMPLITUDE 

In this Appendix we will find the relationship between f and t. We 
start with equation (2.14). Using (2.15) and (2.16), we can write (2.14) 
in terms of the transition operator and Incident field for N particles 

00 

W . I L(L“ 1 V )\ (A. 1) 

8 0 N-l 0 


From (A.l) we can write 

00 

T - l L(L _1 V) N (A. 2) 

N-l 

.From eq. (A. 2), we have 

ty - vy + VL~hrv + vl" 1 vl" 1 v'J' + . .. (a.3) 

Substituting eq. (2.8) for Y, (2.12c) for V and (2.19) for T in eq. (A.3) 
along with the free space dyadic Green’s function for L \ we obtain 


Tf 


T . T 
-.11 *12 

V 


*(x) 0* 

V 


1 (x ) 0 

r 

-11 £l2 

i(x) o" 

E(xO" 



“ 



+ 


1 




T T 

—21 =22 

J 

H 


0 0 

H 


0 0 


_-21 -22 

0 0 

H(x") 


dx:"* 



(A. 4) 
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where % ( x ) ■ iwe ^AU (x) I 


Using eq. (2,18), (A, 4) can be expressed as 



0 

0 


(A. 5) 


Now, we can conclude by comparing both sides of eq. (A, 5) that for dielectric 
scatterer 

«12 * “22 * “21 * ® (A. 6) 

and 



To find the relationship between the transition operator and the 

I 

scattering amplitude, we write out eq, (2.17) explicitly: 


Y s (x) - if 1 ^ - J dx'*G(x,x'*) | dx^t (x",x'"')'F 0 (x"') (A.8) 

Here the free space dyadic Green's function G(x,x") is obtained from 
Felsen and Marcuvitz, Chapter I, Sec. 1.1, it is given by 
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G(x,x') ■ 


■ia)pQ(X +—j) 
k 0 

"Vxl 


7x1 


nn 

-i“*o ( S + 7§> , 
k o J 


4tr | x-x'* | 


(A. 9) 


To obtain 7 (x) in the radiation zone, the far field expression for 
8 “ 

G(x,x ) will be required [Twersky, 1967]: 


G(x,x') * H « 

where 

H - 


-ik n 0‘x" lk 0^ 
0 e 


4it j 3C | 

“iuy 0 (I - 0 0) 
-ik o 0x(I) 


x ■+*> 


ikoPxCp 


-iwe 0 (I- 0 £)J 


(A. 10) 


(A. 11) 


r 

i 


Substituting eqs. (A. 10) and (A. 11) in (A. 8) and using for a incident 
wave Y 0 (x) - Y Q (0) e ^ (A.12) 


and 


i f 0/xO 

tCk^Ojk-i) * r dx'*dx‘*‘'t(x'*,x'*'*) e 

°" °“ (2tt) 3 J ~ ~ 


(A. 13) 


we obtain 


Y s (x) - 2 tt 2 M t(-k 0 g,k () i)f 0 (0) - 


ik 0 >* 


(A. 14) 


where 


¥ - 

s 


r i 

e 

=s 


h 

»s 


(A. 15) 


{ .) 


Using eqs. (A. 15) and (A. 7) and’ substituting H into the eq. (A. 14), we obtain 
for dyadic scattered wave 


e g (x) *’ 2 tt [ (1 - 0 0)i(k () 0,k () i>te 1 (0)] 


ik 0 !— 


(A. 16) 


The definition of scattering amplitude given by [Ishimaru, 1978] 



■ 




g H cx)'Vf (E) (0,l) ± 


i k 0 ln 
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x| •*■ » , f2f ^ 


(A. 17) 


Assuming dyadic incident wave, 6^(0) * (I - i. i), finally we obtain the re- 
quired result 


f (E) (0,i) - 2tt 2 (I " 0 ^-tCk^O.k^Kl -11) 


(A. 18) 


Equation (A. 18) is the same as Lang [1981) obtained with different method. 
From eq. (A. 14) we can also obtain 

ik n 1*1 


h s (x) - 2iTn E (0xl>t (k^, k () i>e i (0) ] 


(A. 19) 


The definition of scattering amplitude for the H field given by 


h (x> ~ I <H) (0,i) - 


ik o'- 


x ■> 00 


(A. 20) 


For a dyadic incident wave e^(0) ■ (I - i i) we have 
| ^ (0,i) * 2tt 2 n (0^x1) •t(k (} 0,k (> i) • (I - 1 i) 


where 


n » * — 

m, 



o f 


(A. 21) 


(A. 22) 


In order to obtain the relationship between two different scattering 
amplitudes, we take the cross product of eq. (A. 21) with respect to £ direction 
we find i 


£ x f^=~2iT 2 n 0 x (Oxl) *£* (I -11) 


but 


0x (0 x I) * - (I -00) 

then eq. (A. 23) becomes 

0 x f <H « n2u 2 (I - 0 0>tr(I - i i) 


(A. 23) 


(A. 24) 


(A. 23) 
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c 

Using eq. (A. 18) In eq, (A. 25) we obtain the following relationship for 


scattering amplitudes for E and H fields 
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APPENDIX B 

EVALUATION OF ZEROTH ORDER SOLUTION 
‘ OF SLOWLY VARYING COEFFICIENTS 


In this section, we will evaluate eq. (3.59), which is given by 

■ o 


where 


4>^(z,z) - - p( 2 )||dz , '|d8|dz > ''y hz^xyCz'")^ 0 ^ (z)| 


rect^(z) 


Performing eq.(]B.l), we obtain 


“<{»^(z) * p(z)M<f>^ (z)re££ ,(z) 
dz a 


where 


Z oo oo 




0 —00 -00 
Now, using eq. (3.51) we lur'e 


M - 


where 


lim 

Z-Ko 


Z 

~ | dz* jdsdz" 

o J 


^“»") X ( “» +) ' 


KH 


^(+»-) ^(+>+) 




X ^^(z'.Z'^B) - ^^X 11 (^,Z-.8,Z^-s)-lV iK(aZ ^ 62 ^ ) 


(B.l) 


(B. 2) 


(B.3) 


(B.4) 


(6.5) 


(B.6) 


a»3 e {+,-} 


The Integrals over ds and dz'"’ will new be evaluated. First define 


g (a,B) - Jdadz''f< a,3) (z‘\z'%a) 

By using eq. (B.6) in eq. (B.7), we obtain 

^ M. t« . Q ( a *P) . 

B ' 2N_ i a m 


where 


g (a » 3) -J dsdz' f '-B)e“ ± ‘ c(az '*~ e?/ ' ) 


The quantity will be related to t^. We note that 

f _ 

-11 ^-t ° 


d *t 


and from eq. (3.86) we have 

Xii(x t ,z'- S ,z~- S ) - jdi t T 11 (-§ t ,-x t: -£ t ,z''-8 > z"- 
We also have from eq. (2.22) 


8) 


t u CXt,*;.z,z') - ^3 JdK £ dK z d]^d^t u fe t , K2 . ! .^, K ') 


" hZ'K + K , z “ O'] 

— t ~t —x. — t z z 

e 


By using eqs. (B.10), (B.ll) and (B.12) in eq! (B.9), evaluating the 
integrals over dx^., ds t »ds,and dz“'" > and using the resulting Dirac delta 
to eliminate the remaining integrals , we obtain 


g (a » 3) - (27r) 3 t 11 (k to ,8tc5k to ,8K)e iK(3 ”° )z " 


Putting eq. (B.13) in eq, (B.8), we have 




m 


(B.7) 

(B.8) 

(B.9) 

(B. 10) 

(B.ll) 


(B.12) 

resulting 

fundtions 

(B.13) 

(B.14) 
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The final expression for M can now be obtained. We put eq, (B.14) into eq. 
(B.5) and perform the integration over z' and also take the limit as z goes to 
infinity. We obtain 

fs" °1 



+ 

in 


(B.15) 


where 



B(2tt) 3 

2^ 




(B.16) 


3 « {+,-> 

Now employing eqs. (2.20), (2.23), (3.27b) and (3.37a), we have 


m 


^ m k 

K - • 


(B.17) 
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APPENDIX C 


EVALUATION OP T AND T 


To find an expression for T we use eq. (3.83) in eq. (3.81). Ws 


1 


A(-d) - - ~-(y ry ) *T*A(-d)e 


KtC -K)d 


(C.l) 


T — 

where (y“ ry g ) is a dyadic and N T -nK/k Q . By taking the inverse of this 
dyadic an expression for T can be obtained. It is: 

_ T _ K<-< )d 

T - =2N T (y Fy g ) -I e 8 (C.2) 

_T _ _i _T _ 

where we have used the identities I * A ■ A and (y Ty ) *(y Ty ) * I# The 

m 8 8 “8 

expression for T can now be obtained by employing eq. (3.84) and eq. (C.2). 


We have 


T T 

ti , + „ -v., - t* “\“1 T “ 2i<d 

£ “ -Cy Ty g )*(y Ty g ) *1 e 


(C. 3) 


+ T _T 

The quantities y ry” and y” Ty” need to be evaluated. From eq. (3.31), 

© o 


we write 


y” - p + p_ h and y » P +1 v“ + ^ +2 h 


, z>-d 


y” « V, v° + h° 
3 1 lg-g 2g - 


{CA) 


where f + ^, ^ lg and asre Slven by eq. (3.28). Using these expressions. 


we find 


y + Ty” ■ cv^v 0 + cv^ h^ + d h°v^ + d h°h° 
3 7 g “+“g -+“ 8 


(C.Sai 


~ V 0 0, *0,0 , ‘ 00,. .0.0 

y Ty • av v„ + av h +b hv+bhh 
g g ”>■— g 


(C.5b) 
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where 



* - (f_ 1 ,r? lg ) i 

■ <*-i’ r W 

(C.6a) 

b - <*-2-*V £ 

■ <’L 2 ’ r V 

(C.6b) 

C - e +1 »rf lg ) ' J 

% 

- < »+l' r ' 1 '2g ) 

(C.6c) 

d - Cf +2 ,ry 2g ) d 


(C«6d) 

Substituting the expressions for 

H^g, i-1,2 in eqs. (C.6) we have 


0 0 ,0 .0. 0 0 
a»nv *2 xh -nh*z xv 
8 - “ “8 

. a - 0 

(C.7a) 

. .0 0 0, 0 0.0 
b - ~r| h *z xv + nv »z xh 
g~ — “g - — — — 

b - 0 

(C.7b) 

00.0 ,0^0 

c * n v . • z xh - nh • z ~ icv 
g~+ — — ~ — -g 

c ■ 0 

(C.7c) 

, , 0 0 0 , 0 0,0 
d « -n h »z xv„ *lr qv , *z xh 
g— “g — h “ — 

3 « o 

(C.7d) 


Now using eqs, (3.24) and (3.25) in the expressions for a-d in eqs. (C.7), 
we have 


a - -n(e g K + < g )/k g 


b - -q(K + K g ) /ICq 


c - n(e gK - K g )/k g 


d - I|(< - Kg) /kg 


(C.8a) 

(C.8b) 

(C.Sc) 

l 

‘ (C.8d) 


Proceeding, we use the orthogonality results of eqs. (C.7) which yield 


+ T „ - 

0 0 

, , , o, 0 

y ry g - 

c V . V 
— +-g 

+ d h h 

_ T 

0 0 

. . 0, 0 

y~ Ty - 
8 

a v v 
g 

+ b jh h 


(C.9a) 

(C.9b) 
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To complete the calculation, the Invirse of y~ Ty^ la required. It 
requlrea that 




<y“ Ty g ) 

- i" 
•8 

(C.10) 

or 


ry g )" X '(*y°_y^ + 

,.0.0. 0 0. .0.0 
hh h ) * vjv + h h 

— — “* ~ 

(C. 11 ) 

If 


, - T „ .-1 10 0 
(y ry g ) - 1 2 ^. 


(C.12) 

then eq.(C.ll) la satisfied. Putting eqB. (G.9a) and (C.12) into eqr). 

» 


(C.2) 

and C.3) 

we have 





±C< -Od 

T ■ T e * 

■ -g 


(C.13) 

where 


0 0 

T » T , h h u + T 
•g gh 

0 0 

VV 
~8 — 

(c.14) 

and 


r - r e 2lKd 

* “g 


(C.15) 

where 


r - r . h°h° + r 

m g gh gv 

0 0 
V V 
— **■ •— 

(C.16) 

The Fresnel transmission and reflection coefficients are given by 



m a 

2< 

- m 2K ^S 



V 

K+K * 

g 

gv Eg * + K g 

(C.17) 

and 

V ' 

K-IC g 

K+K * 

g 

e k - K 

p „ £ £ 

EtC + K 

g 8 

(C.18) 
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CALCULATION OF f^i 0 ,! 0 ) 


In this Appendix, we will show that We start by wrltting 

scattering amplitude that is given by LeVine, ©t al [1982], 


KiV> - 1 $) 


I TAS(v^) (h°h°e£ + v^v 0 ®“)©inc9* 
*' c 


+ (h^h.^ e~ + vP e^) sinc0 + 


(D.l) 


where 1° and 0° are the directions of incident and scattering waves respectively. 
The expressions appearing in eq. (D.l) are defined in chapter V. 

We need in forward direction. Using jL^-0^ and dotting eq. (D.l) by h^ 
from l‘W; right and by v^ from the left, we find: 

f vh “ " 0{[ae*+b + ce*]sinc0”’+ [ae” + bj:e“]sinc0 + } (D.2) 


where 


(fe)“ 5 o 


(D.3) 


, 0 . 0 . ,,0 . 0 . 

a " ) (h 

. 0 
V M + 

- X e - 

0 , - 
c * v *h, 



and h^, v® are the unit horizontal and vertical vectors of the incident wave. 

Next by expanding the unit vectors in cartesian coordinates, we determine 
that f ^ has the following dependence on $ * $ - 
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f y k * h(sin$)cos^ (D.4) 

whet a h(c) la a single valued function of 5. Using «q. (D.4), and ass using 

$ la uniformily distributed, we have: 

TT TT * . TT 

? hv - £ I V* ’ If f - u I Krt&'o&f - 0 <»■*> 

“TT -W -TT 


The last integral is zero because of the even and odd character of 
cos <p and sin $j thus f^^O* We also find that **^"0 using the same method. 
Using a similar procedure one can show that 

+ 

K 


(D.6) 
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